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A multicast  group  is  a set  of  nodes  that  are  the  common  destinations  of  the  same 
group  of  messages.  The  source  or  sources  may  be  nodes  within  the  multicast  group  or 
may  be  other  nodes  in  the  network.  Each  message  is  addressed  to  all  the  elements  of 
the  group.  Some  applications  require  that  the  multicast  messages  be  received  in  the 
same  relative  order  by  the  nodes  in  the  multicast  groups.  Our  objective  is  to  provide 
the  multiple  group  ordering  property  for  some  known  topologies  (linear,  ring,  mesh, 
hypercubes)  by  constructing  multicast  trees.  Two  different  graph  models  can  be  used 
to  guarantee  the  message  ordering  properties.  In  the  first  one,  a single  global  tree  is 
built  for  all  the  nodes  in  overlapping  groups,  while  in  the  second  one,  single  trees  are 


built  only  for  overlapping  nodes,  and  on  top  of  these,  separate  trees  are  built  for  each 


multicast  group.  The  goal  in  both  methods  is  to  find  multicast  trees  which  minimize 
time  and  traffic. 


IX 


CHAPTER  1 
INTRODUCTION 


A muIticcLst  group  is  a set  of  nodes  that  are  the  common  destinations  of  the  same 
group  of  messages.  The  source  or  sources  may  be  nodes  within  the  multicast  group, 
or  may  be  other  nodes  in  the  network.  Each  message  is  addressed  to  all  the  nodes  of 
the  group.  Broadcasting  is  a special  case  of  multicasting,  where  the  nodes  that  make 
the  group  are  all  the  nodes  in  the  system. 

There  are  many  applications  in  which  it  is  necessary  that  the  nodes  in  a multicast 
group  receive  two  or  more  messages  in  the  same  relative  order.  Suppose,  for  exam- 
ple, that  nodes  A^i  and  N2  receive,  from  the  same  source  or  from  different  sources, 
messages  mi  and  m2.  Nodes  A^i  and  N2  should  order  the  messages  mi  and  m2  in 
the  same  way.  That  is,  the  multicast  protocol  must  guarantee  that  both  nodes  agree 
upon  the  order  in  which  the  messages  will  be  processed.  Consider  for  example  a bank 
with  two  main  computers.  Each  of  these  maintains  a copy  of  the  banking  database, 
and  they  will  process  transactions  from  the  branch  offices.  In  this  case,  the  two  main 
computers  may  constitute  a multicast  group,  and  each  branch  office  is  a potential 
source.  Now  consider  a deposit  and  a withdrawal  to  the  same  account.  If  the  with- 
drawal is  done  first,  an  overdraft  may  occur  and  a penalty  charged.  If  the  deposit  is 
done  first  there  may  be  no  penalty.  Thus  a different  account  balance  may  result  if 
the  two  transactions  are  not  executed  in  the  same  relative  order. 
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Three  different  ordering  properties  can  be  required  by  a multicast  protocol: 

• single  source  ordering  - in  this  case  messages  will  be  processed  in  the  same 
relative  order  only  if  they  are  originated  by  a single  source,  and  are  destined  to 
the  same  multicast  group. 

• multiple  source  ordering  - in  this  case  messages  will  be  processed  in  the  same 
relative  order  even  when  they  originate  from  different  sources,  but  are  addressed 
to  the  same  multicast  group. 

• multiple  group  ordering  - this  provides  the  general  condition,  in  which  we  can 
have  multiple  sources,  and  the  messages  can  be  processed  in  the  same  relative 
order  even  when  they  are  addressed  to  different  but  overlapping  groups. 

The  banking  example  given  above  must  satisfy  property  b).  Now  consider  adding 
a second  multicast  group  to  distribute  new  software  releases  or  new  system  tables 
(like  defining  overdraft  penalty  charges).  This  second  multicast  group  consists  of  the 
two  main  computers  and  other  development  machines.  It  may  be  required  that  the 
messages  still  be  processed  in  the  same  relative  order  by  the  intersection  of  the  two 
multicast  groups.  In  this  case  multiple  group  ordering  must  be  provided. 

It  is  easy  to  satisfy  single  source  ordering.  The  source  can  cissign  a sequence 
number  to  every  message,  providing  a total  ordering  of  the  messages.  It  is  harder  to 
deal  with  the  other  two  properties.  In  the  second  chapter,  existing  algorithms  to  solve 
this  problem  are  reviewed.  The  third  chapter  contains  a description  of  our  approach 
to  attack  the  problem.  Chapters  four,  five,  six  and  seven  describe  the  results  we 
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have  obtained  for  linear,  ring,  mesh  and  hypercube  networks,  respectively.  The  last 
chapter  concludes  our  work,  and  points  out  directions  for  future  research  on  this  area. 


CHAPTER  2 
RELATED  WORK 


Kaaishoek  et  al.  [26]  have  proposed  an  atomic,  reliable  broadcast  protocol  for 
a broadccist  medium  based  on  a central  site  that  assumes  the  role  of  a central  se- 
quencer. When  a process  wants  to  broadcast  a message,  it  sends  a point-to-point 
message  to  the  central  sequencer.  The  sequencer  assigns  the  next  sequence  number 
to  the  message  and  broadcasts  the  message.  In  this  way,  the  sequencer  provides  a 
global  order  for  every  broadcast  message.  All  broadcasts  are  issued  by  the  sequencer. 
The  protocol  handles  missing  messages  in  a simple  manner,  and  with  low  overhead. 
The  sequencer  maintains  a history  buffer  to  keep  track  of  the  messages  that  have  been 
received  by  every  node.  The  nodes  communicate  this  information  to  the  sequencer 
in  a header  field  every  time  they  send  a point-to-point  message  to  the  sequencer  for 
a broadcast  request.  This  way  the  sequencer  maintains  a table,  indexed  by  node 
number,  containing  the  information  that  node  i has  received  at  least  the  messages 
from  0 to  T{.  The  sequencer  can  compute  the  lowest  value  in  the  table  and  discard 
all  messages  up  to  that  value,  since  it  knows  that  all  nodes  have  received  such  mes- 
sages. The  protocol  requires  two  messages  per  broadcast  in  most  cases.  A major 
disadvantage  of  this  method  (and  of  any  centralized  method),  however,  is  that  the 
central  site  could  become  a bottleneck  in  the  system. 
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Another  centralized  approach  has  been  proposed  by  Chang  and  Maxemchuk  [4]. 
As  in  the  previous  protocol,  the  approach  is  based  on  a central  site  that  provides  the 
global  ordering  of  messages.  The  central  site  is  identified  with  a token.  However,  in 
this  protocol  the  role  of  central  site,  or  token  site,  is  rotated  among  all  the  nodes 
in  the  network.  The  source  sends  the  broadccist  message  to  all  the  nodes,  and  the 
token  site  acknowledges  the  message.  The  positive  acknowledge  contains  a sequence 
number  to  provide  the  message  ordering.  On  each  acknowledgment  the  token  is 
transferred  to  the  next  site.  A site  accepts  the  token  only  if  it  has  received  all  the 
messages  so  far.  This  provides  fault-tolerance  to  the  protocol.  Tokens  are  transferred 
and  accepted  in  acknowledgments,  so  they  in  general  do  not  generate  additional 
messages  if  there  are  messages  to  be  acknowledged.  Depending  on  system  utilization, 
however,  they  may  generate  three  messages  per  broadcast:  the  broadcast  message, 
the  positive  acknowledge  sent  by  the  token  site,  and  the  transfer  of  the  token  after 
each  acknowledgment.  Fault-tolerance  is  obtained  by  transferring  the  token.  If  L 
nodes  accept  the  token  after  a message  hcis  been  acknowledged,  then  the  protocol 
can  tolerate  L processor  failures.  Long  delays  are  introduced,  however,  before  the 
messages  can  be  processed. 

Some  distributed  approaches  have  also  been  proposed.  The  first  widely  used 
method  to  establish  order  among  events  in  a distributed  system  was  proposed  by 
Leslie  Lamport  [27].  In  this  approach,  timestamps  are  assigned  to  every  message  at 
the  source,  and  the  messages  are  delivered  in  timestamp  order.  Using  logical  clocks 
in  the  sites  a total  ordering  of  the  messages  can  be  obtained.  The  major  problem 
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with  this  kind  of  approach  is  that  every  process  must  know  about  the  other  process’s 
operations.  To  guarantee  the  multiple  source  and  multiple  group  ordering  properties, 
when  a node  receives  a message  with  its  timestamp,  before  processing  that  message 
the  node  must  check  all  nodes  in  the  network  (all  of  them  are  potential  sources  to 
the  multiccLst  group)  for  other  messages  not  yet  received  with  smaller  timestamps 
addressed  to  the  same  multicaist  group  or  to  an  overlapping  group.  Only  when  the 
node  is  certain  that  it  has  received  the  message  with  the  smaller  timestamp  can  it 
proceed  to  process  the  message. 

Birman  and  Joseph  [3]  have  proposed  a two-phase  protocol  to  provide  total  order- 
ing of  messages.  Each  site  maintains  a priority  queue  cissociated  with  each  process. 
When  a node  receives  a message,  it  assigns  a unique  priority  number  to  it  higher 
than  any  other  given  to  that  process  before,  marks  the  message  as  undeliverable, 
and  sends  this  priority  number  to  the  sender  of  the  message.  The  sender  computes 
the  maximum  priority  number  from  all  the  priorities  received  from  all  active  desti- 
nations, and  sends  this  value  back  to  the  destinations.  These  change  their  initial 
priority  number  to  the  new  one,  tag  the  message  as  deliverable,  and  re-arrange  their 
priority  queues.  A message  at  the  front  of  the  queue  whose  tag  is  deliverable  is  then 
delivered.  Their  protocol  is  mainly  concerned  with  fault  tolerance,  and  hence  they 
may  have  a large  message  overhead. 

A different  approach  is  followed  by  Garcia-Molina  and  Spauster  [13].  For  a given 
set  of  multicast  groups  with  overlapping  nodes,  a single  tree  is  generated,  called  a 
propagation  graph.  The  graph  indicates  the  path  that  each  message  must  follow 
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to  get  to  all  the  destination  nodes.  The  sender  sends  the  message  to  a primary 
destination  of  the  multiccist  group,  which  is  the  smallest  common  ancestor  to  all  the 
nodes  in  the  group  in  the  propagation  graph.  Instead  of  ordering  the  messages  at  a 
single  central  site,  they  are  ordered  by  the  collection  of  nodes  in  the  tree  by  merging 
messages  destined  for  different  groups.  Nodes  that  are  in  the  intersections  of  the 
multicast  groups  are  selected  to  be  the  intermediary  nodes  in  the  graph. 

One  technique  that  has  been  applied  to  provide  multicaisting  in  computer  networks 
is  single-tree  and  multiple-tree  forwarding.  In  this  technique,  single  or  multiple  span- 
ning trees  are  constructed  that  span  the  nodes  in  the  multiccist  group.  Source-b2ise 
forwarding  may  be  used  to  construct  such  a tree  for  each  node  in  the  system.  These 
trees  can  minimized  both  delay  and  bandwidth.  The  algorithms  proposed,  however, 
generally  do  not  consider  multiple  groups. [11] 

Minimum  spanning  trees  and  shortest  path  trees  have  been  proposed  to  build 
trees  for  broadcasting  and  multicasting.  Minimum  spanning  trees  [9]  minimize  the 
total  branch  cost  for  a broadca.st.  For  multicasting,  minimum  Steiner  trees  must 
be  constructed  to  minimize  the  total  cost  of  a multicast.  However,  finding  minimal 
Steiner  trees  has  been  shown  to  be  an  NP-Complete  problem  [25].  In  any  case,  these 
trees  in  general  do  not  minimize  message  delay.  Shortest  path  trees  [37]  minimize 
the  distance  from  the  root  to  any  other  node  in  the  tree.  Wall  [37]  discusses  several 
ways  of  computing  ”a  center”  of  a network  and  constructing  its  shortest  path  tree, 
so  that  the  trees  constructed  have  good  delay  properties. 
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Lan  et  al.  [28]  considered  the  problem  of  constructing  shortest  path  trees  and  at 
the  same  time  minimizing  traffic  for  multicasting  in  hypercube  networks.  Esfahanian 
and  Choi  [6]  have  shown  that  building  optimal  multicast  trees  is  an  NP-Complete 
problem  for  hypercubes.  They  developed  a heuristic  algorithm,  taking  cis  the  root  of 
the  tree  the  source  of  the  multicast  message.  This  approach  leads  to  a multicast  tree 
for  each  different  source,  which  is  inappropriate  for  message  ordering. 

The  Steiner  tree  problem  is  closely  related  to  the  Optimal  Multicast  Tree  problem 
(OMT).  An  OMT  for  a multicast  group  is  defined  as  a tree  that  minimizes  the  max- 
imum distance  from  the  root  to  the  farthest  element  from  the  root  in  the  multicast 
group  that  has  minimum  traffic,  that  is,  minimum  number  of  links.  A minimum 
Steiner  tree  minimizes  only  the  traffic  parameter.  In  general,  optimal  multicast  trees 
do  not  minimize  traffic,  i.e.,  they  are  not  minimum  Steiner  trees.  Garey  and  John- 
son [15]  have  shown  that  finding  minimum  Steiner  trees  for  graphs  with  rectilinear 
distances  is  an  NP-Complete  problem.  This  has  lead  us  to  conjecture  that  building 
an  OMT  for  mesh  networks,  which  use  rectilinear  links,  is  also  NP-Complete.  The 
problem  of  building  an  OMT  in  a grid  for  a multicast  group  D given  the  root  is 
equivalent  to  the  problem  of  finding  a rectilinear  directed  Steiner  tree  for  D in  a 
grid  in  which  all  the  links  point  away  from  the  root  node.  A special  case  of  this 
problem  (in  which  all  the  nodes  lie  in  the  first  quadrant  of  and  the  directed  tree 
is  rooted  at  the  origin)  has  been  known  in  the  literature  as  the  Rectilinear  Steiner 
Arborescence  problem  (RSA).  Rao  et.  al  [32]  have  proposed  a heuristic  algorithm 
for  this  problem  that  has  time  complexity  of  0(|  D | /o^r  j £)  |).  They  proved  that 
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the  traffic  generated  by  the  heuristic  has  an  upper  bound  of  twice  of  that  generated 
by  any  optimal  algorithm.  They  have  proposed  an  extension  of  this  algorithm  that 
solves  the  general  case  in  which  the  nodes  in  D lie  anywhere  in  the  plane  that  runs 
with  0(1  D p log  \ D ])  time  complexity.  We  propose  in  this  dissertation  a similar 
heuristic  with  same  upper  bound  that  also  has  time  complexity  0(|  D \ log  | D |) 
but  solves  this  general  case  in  which  the  nodes  are  allowed  to  be  anywhere  in  E^. 


CHAPTER  3 

OPTIMAL  MULTICAST  TREES 


The  objective  of  our  work  is  to  present  a technique  to  guarantee  the  message 
ordering  properties  for  several  known  topologies:  linear,  ring,  mesh  and  hypercubes. 
Multicast  messages  will  be  ordered  by  a collection  of  nodes  structured  in  a message 
propagation  graph.  The  graph  indicates  the  flow  of  the  messages  through  the  network. 

An  approach  that  hcis  been  described  to  handle  multiccisting  is  the  single-tree 
forwarding  technique.  In  this  technique  a multicast  tree  is  constructed  that  spans 
the  nodes  in  the  multicast  group.  Each  node  sends  message  copies  along  the  branches 
of  the  tree.  Given  a single  multicast  group,  a multicast  tree  can  guarantee  that  all  the 
nodes  in  the  group  receive  the  messages  in  the  same  relative  order.  The  root  of  the 
tree  decides  the  message  ordering.  With  multiple  groups,  single  trees  are  needed  that 
span  the  nodes  in  each  multicast  group.  To  guarantee  the  multiple  group  ordering 
property,  a single  tree  (or  subtree)  can  be  constructed  for  overlapping  nodes  in  the 
multicast  groups.  Two  different  graph  models  could  be  used  to  guarantee  the  message 
ordering  properties.  In  the  first  one,  a single  global  tree  is  built  for  all  the  nodes  in  the 
overlapping  groups,  while  in  the  second  one,  single  trees  are  built  only  for  overlapping 
nodes,  and  on  top  of  this  subtrees,  separate  trees  are  built  for  each  multicast  group. 
In  any  of  these  graphs,  the  source  sends  a multicast  message  to  the  closest  ancestor 
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of  all  the  nodes  in  the  multicast  group.  This  common  ancestor  is  called  the  primary 
destination  of  the  multicast  group. 

It  is  assumed  that  any  site  can  be  the  source  of  messages  to  a particular  destination 
group.  When  a node  wants  to  send  a message  to  a particular  group,  a path  from  the 
source  to  the  primary  destination  must  be  determined.  In  the  process  of  constructing 
multicast  trees,  the  source  of  the  messages  is  not  taken  into  consideration,  since  it  is 
assumed  that  the  sources  to  the  multicast  groups  are  uniformly  distributed. 

Two  parameters  must  be  considered  to  construct  a multicast  tree:  time  and  traffic. 
The  parameter  time  is  a measure  of  the  time  required  to  send  a message  from  the 
primary  destination  of  a multicast  group  to  the  last  node  in  the  group  that  receives 
the  message.  The  time  to  send  a message  from  one  node  to  another  is  given  by  the 
distance  between  the  two  nodes.  The  cost  of  the  multicast  is  associated  with  the  time 
required  to  complete  the  multicast,  and  is  determined  by  the  distance  between  the 
primary  destination  and  the  node  in  the  group  farthest  from  the  primary  destination. 
A message  can  be  forwarded  in  parallel  through  disjoint  paths  in  the  tree.  The 
traffic  parameter  is  given  by  the  number  of  links  needed  to  multicast  the  message 
from  the  primary  destination  to  every  element  of  the  group  through  the  branches 
of  the  multicast  tree.  The  messages  should  arrive  as  fast  as  possible  to  all  their 
destinations,  using  as  few  links  as  possible  to  achieve  that.  These  two  parameters 
are  not  independent.  For  example,  consider  a mesh  network  in  which  the  nodes 
(2, 1),  (3, 1),  (3, 2),  (3, 3),  (2, 3)  constitute  a multicast  group  and  the  root  of  the  tree 
is  given  by  the  node  (1,1).  The  multicast  tree  that  minimizes  the  traffic  generated 
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is  given  by 

(1,1)  ^ (2,1)  ^ (3,1)  ^ (3,2)  ^ (3,3)  (2,3) 

This  tree  generates  a traffic  of  5 units.  However,  the  cost  (or  time)  is  also  5,  while 
the  tree  given  by 

(1.1) 

(ii) 

/\ 

(2.2)  (3,1) 

I I 

(2.3)  (3,2) 

I 

(3,3) 

generates  a traffic  of  6,  but  its  cost  is  only  4 units. 

An  optimal  multicast  tree  (OMT)  is  a multicast  tree  with  a minimum  cost  that 
hcLS  the  lowest  traffic. 

3.1  Global  Single  Tree 

Suppose  we  partition  the  groups  so  that  no  group  in  a partition  intersects  a group 
in  another  partition,  and  we  form  as  many  partitions  as  we  can  with  this  property.  A 
single  tree  can  be  constructed  that  spans  the  nodes  in  each  partition.  For  example, 
consider  the  following  multicast  groups:  A={a,b,c,d},  B={d,e,f},  C={a,b,f},  where 
the  capital  letters  A,B,C  represent  the  multicast  groups,  and  the  lower  case  letters 
represent  the  nodes.  Suppose  that  the  connections  of  the  network  are  given  by  the 
following  graph: 
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b 


Then  an  optimal  multicast  tree  is  given  by: 


b 


Node  b is  the  primary  destination  for  groups  A and  C,  and  node  d is  the  primary 
destination  for  group  B.  Two  different  messages  sent  to  groups  A and  C are  ordered 
by  node  b. 

Lemma  3.1  The  multiple  group  ordering  property  is  satisfied  under  the  global  single 
tree  model. 

Proof: 

We  need  to  show  that  any  two  nodes  a and  b belonging  to  two  or  more  groups  will 
receive  the  messages  addressed  to  these  groups  in  the  same  relative  order.  Since 
all  of  these  groups  contain  nodes  a and  b,  it  must  be  the  case  that  their  primary 
destinations  lie  on  or  above  a and  b in  the  multicast  tree.  Let  P be  the  lowest  one 
of  these  primary  destinations,  i.e.,  the  closest  common  ancestor  to  a and  b.  Every 
message  addressed  to  any  of  the  groups  that  contain  nodes  a and  b ha.s  to  go  through 
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node  P before  reaching  nodes  a and  h.  The  order  of  the  messages  for  the  nodes  in 
the  groups  that  contain  both  a and  h can  thus  be  decided  by  node  P,  and  nodes  a 
and  b will  agree  in  that  order.  □ 

Essentially,  two  steps  are  required  to  construct  the  multicast  tree  for  a partition. 
First,  an  appropriate  root  must  be  found.  Once  an  optimal  root  is  computed,  we  must 
determine  the  best  routing  strategy  to  cover  every  node  in  the  multicast  groups,  and 
determine  the  primary  destinations. 

Consider  first  the  problem  of  finding  an  optimal  root.  We  assume  that  sources 
are  uniformly  distributed  among  the  nodes  of  the  network.  We  also  cissume  that 
the  probability  that  an  arbitrary  message  is  addressed  to  a particular  destination 
group  is  the  same  for  every  group.  We  want  to  select  the  root  such  that  the  sum  of 
all  distances  from  the  root  to  the  elements  in  each  group  farthest  from  the  root  is 
minimized.  The  problem  can  be  formulated  as  follows. 

Let  N be  the  number  of  nodes  in  the  network.  The  set  of  nodes  will  be  represented 
by  the  set  {0, 1, . . . , — 1}.  Let  d(x,  y)  be  the  length  of  a shortest  distance  between 

nodes  x and  y.  Let  be  the  jth  node  of  the  ith  group  in  a multicast  partition. 

The  optimal  root  of  the  OMT  is  defined  as  the  node  that  produces  the  minimal 
cost.  The  cost  for  a root  node  i is  defined  by: 

P 

c,  = , ^max  d{pdi,Tijk) 

j — I i S K < TTlj 

where  p is  the  number  of  multicast  groups,  pd,  is  the  primary  destination  for  group 
i,  and  is  the  number  of  nodes  in  the  jth  multicast  group.  The  root  node  with 
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minimal  cost  Cm,  is  thus 


min  c, 
0 <i < N -1 


A problem  here  is  that  the  primary  destinations  are  not  known  before  the  tree  is 
built,  and  even  before  the  root  has  been  found.  In  some  ca^es  it  may  be  necessary  to 
find  the  root  that  minimizes  the  sum  of  the  distances  from  the  root  to  the  farthest 
element  in  each  multicast  group,  and  not  necessarily  the  sum  of  the  corresponding 
distances  from  the  primary  destinations.  In  this  case  we  would  be  minimizing  the 
maximum  potential  cost  of  the  multica.st  tree,  and  not  the  actual  cost  of  the  tree. 
The  cost  of  the  tree  is  then  given  by 


P 

Ci  — ^ max  d(i,rijk) 
j — ll  < k < rrij 

A naive  way  of  computing  the  cost  of  every  node  requires  time  complexity  0{Nm), 
where  m = rrij,  assuming  that  the  distances  between  each  pair  of  nodes  can  be 
computed  in  constant  time. 

Once  we  have  the  root,  the  rest  of  the  tree  must  be  determined.  Routes  must 
be  found  from  the  root  to  each  node  in  the  multicast  groups.  Determining  the  best 
routing  for  linear  and  ring  networks  is  simple  and  can  be  done  in  time  proportional 
to  N . For  mesh  and  hypercubes  networks  the  problem  is  much  harder.  Shortest-path 
trees  will  be  considered  for  these  networks.  A shortest-path  tree  is  a tree  that  spans 
all  the  elements  of  the  destination  lists,  with  the  property  that  the  path  from  the 
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root  to  every  node  in  the  tree  is  a shortest  path  in  the  original  network.  The  routing 
problem  can  be  formulated  as  follows. 

Let  G be  the  graph  determined  by  the  network  topology,  and  let  dG{x,y)  be  the 
length  of  a shortest  path  in  graph  G from  node  x to  node  y.  We  want  to  construct  a 
tree  T of  root  node  R with  the  following  properties: 

• V{T)  C V{G) 

• the  sum  of  the  maximum  distances  from  the  root  node  R to  the  element  in  each 
multicast  group  farthest  from  R is  minimized 

• (It{x,  R)  = (1g{x,  R)Wx^T 

• no  tree  exists  satisfying  the  three  properties  above  that  has  fewer  nodes  than 
this  one. 

3.2  Multiple  Trees 

The  multiple  group  ordering  property  can  be  satisfied  without  constructing  a 
single  global  tree,  as  long  as  overlapping  nodes  belonging  to  two  or  more  groups 
belong  to  common  subtrees  in  the  trees  for  each  of  these  multicast  groups.  Given  the 
multicast  groups,  define  graph  G'  = {V\E')  such  that  V'  is  the  set  of  nodes  in  the 
groups,  and  for  Vi,V2  G G , the  edge  (i>i,U2)  ^ E iff  G f)  gj  for  any  two 

groups  gi  and  gj.  That  is,  two  nodes  are  connected  in  G'  if  and  only  if  they  both 
belong  to  the  intersection  of  at  least  two  groups.  If  (ui,U2)  G E'  then  vi  and  are 
overlapping  nodes,  so  they  must  receive  messages  addressed  to  groups  i and  j in  the 
same  relative  order.  They  will  belong  to  a common  subtree  of  the  multicast  trees  for 
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these  groups.  In  fact,  all  connected  components  of  G will  be  in  common  subtrees. 
For  example,  consider  the  following  groups:  A={a,b,c,g},  B={b,c,d),  C={c,d,e,f}, 
and  D={e,f}.  Then  the  graph  G'  is  given  by: 


a 

b 


A multiccLst  tree  T,  has  to  be  constructed  for  every  group  i.  The  next  step  is 
to  build  subtrees  that  span  the  nodes  in  each  connected  component  of  G" , and  then 
connect  these  subtrees  such  that  the  following  condition  is  satisfied:  a subtree  for  a 
connected  component  would  be  a subtree  of  T,  if  the  connected  component  contains 
at  least  two  elements  of  group  i.  In  the  previous  example,  the  subtree  with  the 
connected  component  {b,c,d}  must  be  a subtree  of  the  multicast  tree  for  groups  A, 
B and  C,  and  the  one  with  {e,f}  must  be  a subtree  for  groups  C and  D. 

Once  the  subtrees  have  been  built,  a multicast  tree  T",  is  built  for  each  group  i in 
the  following  way: 

1.  First  we  need  to  determine  which  subtrees  Sj  need  to  be  part  of  T{.  Every  leaf 
in  Sj  must  belong  to  group  i.  The  nodes  and  links  in  the  subtree  that  do  not 
have  descendants  in  group  i are  not  included  in  the  subtree  for  T,.  Let  A be  the 
set  of  nodes  to  be  connected  in  step  3.  Add  the  root  of  Sj  to  set  A,  and  assign 
a depth  to  this  node  equal  to  the  depth  of  Sj.  In  case  this  node  was  assigned 
a previous  depth  a:  > 0,  then  assign  to  it  maxja:,  depth  ofSj}. 
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2.  Consider  the  nodes  j in  group  i that  did  not  belong  to  any  connected  compo- 
nent. If  node  j can  be  covered  by  any  of  the  subtrees  found  in  the  previous  step 
without  increasing  the  cost  of  the  subtree,  then  considere  this  node  as  covered 
by  the  subtree  (this  may  affect  the  traffic  generated  by  the  subtree,  but  not  its 
cost).  This  node  is  inserted  below  the  closest  node  to  it  already  in  the  subtree 
that  keeps  the  subtree  a shortest-path  tree.  The  longest  possible  path  already 
in  the  tree  is  used  as  a shortest-path  from  the  root  to  node  j.  This  helps  in  re- 
ducing the  traffic.  The  other  nodes  that  cannot  be  inserted  without  increasing 
the  cost  of  the  subtree  are  added  to  A and  a depth  of  0 is  assigned  to  them. 

3.  Connect  the  nodes  in  set  A. 

Consider  the  multicast  groups  given  above,  together  with  the  following  network: 


The  subtrees  generated  for  the  two  connected  components  may  be  as  follows: 


Now  we  need  to  determine  which  of  these  subtrees  (completely  or  partially)  must 
belong  to  the  multicast  tree  for  each  group.  These  subtrees  must  be  connected 
together,  and  to  the  extra  nodes  in  the  group  that  did  not  belong  to  any  connected 


component. 
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For  group  A,  for  example,  only  the  first  subtree  must  be  part  of  its  multicast  tree. 
Nodes  a and  g must  be  connected  to  this  subtree.  A final  multicast  tree  for  this 
group  could  be 

a 

/ \ 

g d 


b c 


For  groups  B,C  and  D the  following  trees  could  be  obtained: 

C D 

Ar 

C * 

The  primary  destination  for  group  A is  node  a,  for  group  B is  d,  and  for  groups 
C and  D is  node  e. 

Lemma  3.2  The  multiple  group  ordering  property  is  satisfied  under  the  multiple  trees 
model. 

Proof: 

We  need  to  show  that  two  nodes  a and  b belonging  to  two  or  more  groups  will  receive 
the  messages  addressed  to  these  groups  in  the  same  relative  order.  Suppose  that 
nodes  a and  6 belong  together  to  the  groups  in  set  C = | a,6  G ^,).  Nodes  a 

and  b will  be  connected  by  a direct  link  in  graph  G' . Therefore,  they  will  belong 
to  a common  subtree,  and  this  subtree  will  be  part  of  the  multicast  tree  for  each  of 
the  groups  in  C.  Let  P be  the  root  of  this  subtree.  Every  message  addressed  to  the 
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groups  in  C must  go  through  node  P before  reaching  nodes  a and  b.  Thus  P can 
provide  the  ordering  of  the  messages  for  the  nodes  belonging  to  the  groups  in  C.  □ 

Algorithms  to  construct  single  trees  given  a set  of  nodes  are  given  in  this  disserta- 
tion for  linear,  ring,  mesh  and  hypercube  networks.  These  algorithms  will  be  used  to 
generate  the  subtrees  for  each  connected  component  of  G when  we  want  to  construct 
multiple  trees.  They  are  used  as  well  to  connect  the  necessary  subtrees  together  and 
the  rest  of  the  nodes  of  each  group  that  did  not  belong  to  any  connected  component. 

In  the  single  global  tree  case,  an  optimal  root  for  the  multicast  tree  is  determined 
given  the  members  of  the  multicast  groups.  Here,  an  optimal  root  must  be  found 
given  a set  of  subtrees  and  a set  of  single  nodes.  The  subtrees  will  be  connected  via 
their  roots.  To  compute  the  optimal  root  in  this  case,  the  previous  algorithms  would 
be  modified  to  take  into  consideration  the  depth  of  the  subtrees  to  be  connected.  Each 
root  of  a subtree  is  cissigned  a depth  equal  to  the  depth  of  the  subtree  it  represents, 
and  the  single  nodes  that  did  not  belong  to  any  intersections  are  assigned  a depth 
of  0.  The  whole  idea  is  again  to  minimize  the  cost  of  the  multicast  for  each  group. 
The  structure  of  the  subtrees  generated  by  the  connected  components  would  remain 
intact  to  guarantee  the  multiple  group  ordering  property.  Once  the  optimal  root  is 
found,  the  routing  from  this  root  to  the  nodes  connected  is  determined.  The  links 
used  in  these  routes  may  include  some  links  already  used  in  the  subtrees.  Thus  it 
may  happen  that  a message  has  to  go  through  one  or  more  nodes  more  than  once 
until  the  multicast  is  completed.  This  fact  tends  to  increase  the  amount  of  traffic 
generated  as  compared  to  the  single  global  tree  solution,  but  the  cost  of  the  multicast 
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would  be  reduced,  and  much  more  variability  in  the  primary  destinations  would  be 
obtained. 

In  the  following  chapters,  the  two  graph  models  to  construct  multicast  trees  are 
discussed  in  more  detail  for  each  of  the  network  topologies  considered,  and  their 
performance  evaluation  compared  in  terms  of  the  costs  of  the  multica.st,  the  amount 
of  traffic  generated,  and  in  terms  of  the  percentage  of  multica.st  groups  that  have  the 
same  primary  destination. 


CHAPTER  4 

OMT  FOR  LINEAR  NETWORKS 
4.1  Single  Group 

Let  the  multicast  group  be  composed  by  the  nodes  ni , . . . , n^.  Let  / = mini<,<„  n,, 
and  h = maxi<,<m  n;.  Assuming  that  the  multicast  group  is  sorted  in  increasing  or- 
der, then  I = rii,  and  h = n^.  The  optimal  root  r is  computed  by  r = [(/  -|-  h)/2\ 

Clearly,  r is  the  midpoint  of  the  interval  from  ni  to  n/,,  and  = mino<,<A/_i  c,- 
The  optimal  root  can  be  found  in  0(1).  In  a linear  network,  constructing  the 
multicast  tree  once  we  have  the  optimal  root  is  trivial.  There  is  only  one  path  from 
each  node  to  every  other  node,  so  there  is  no  room  to  minimize  traffic.  This  unique 
path  has  to  be  taken  in  the  multicast  tree,  and  only  one  message  is  sent  by  every 
node  in  the  tree  to  each  child. ^ 

Example  J.l  : 

Let  N = 8,  and  the  multicast  group  be  composed  of  the  nodes  3,6,  and  7.  Then  the 
optimal  root  is  node  5,  and  the  optimal  multicast  tree  is; 

‘Actually,  every  node  except  the  root  has  only  one  child  in  the  tree. 
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5 

\ 

4 6 

3 7 

If  a node  has  to  multicast  a message  to  this  group,  it  sends  the  message  to  node 
5 through  a shortest  path,  and  the  message  gets  propagated  down  the  tree  to  the 
nodes  in  the  destination  list.  In  this  case,  we  say  that  the  cost  of  the  multicast  is  2, 
which  is  the  cost  of  the  root.  □ 

4.2  Multiple  Groups 

Let  p be  the  number  of  multicast  groups.  Let  /,  = mini<j<m,  and  /i,-  = 
maxi<j<„i,  n,j.  Compute,  VI  < i < p,  r,  = (/,  + /i,)/2.  Each  of  this  p values  represent 
the  midpoint  of  each  group.  The  following  algorithm  finds  the  optimal  root.  Note 
that  the  values  of  the  r'-s  are  either  integers  or  lie  exactly  between  two  integers. 


Algorithm  i./  Linear- Multiple  Groups 
begin 

for  i=l  to  p do  r,  = {L  + /i,)/2; 
sort  the  r'-s] 
if  p is  odd  then 
m = (p+  l)/2; 

if  rm  is  an  integer  then  root  = r^ 
else 

let  n/=number  of  r,’s  strictly  less  than  r„,; 
let  n/?,=number  of  r,’s  strictly  higher  than  I'm', 
if  nl  < nh  then  root  = 
else  root  = 


else 
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m = p/2; 

if  Tm  is  an  integer  then  root  = 
else  if  r„i  < r^+i  then  root  = [r^] 
else 

let  nl  and  nh  be  as  before; 
if  nl  < nh  then  root  = 
else  root  = [r^J 

end; 

Proof  of  correctness: 

There  are  essentially  4 different  cases  we  need  to  consider  to  prove  that  the  optimal 
root  is  computed  by  this  algorithm.  We  will  show  that  for  all  of  them  the  computed 
value  is  correct. 

• Case  1:  In  this  case,  p is  odd  and  is  an  integer,  is  exactly  the  midpoint  in 
the  ordered  list  ri, . . . , Tp.  Consider  the  cost  of  a node  x,  such  that  x + 2 = r„j, 
2 > 0.  The  number  of  groups  with  midpoint  r,  < x is  less  than  50%  of  the 
number  of  groups  (less  than  p/2).  For  this  groups  the  cost  of  choosing  x as  the 
root  will  decrecLse  by  at  most  z units,  as  compared  to  the  cost  of  r^,  while  for 
more  than  50%  of  the  groups  the  cost  will  increase  by  z.  The  same  argument 
applies  for  the  case  in  which  x > r^-  Thus  Cr„,  is  necessarily  optimal. 

• Case  2:  p is  odd  and  is  not  an  integer.  For  this  case,  we  have  computed  the 
value  (r„,)  for  which  the  cost  is  optimized,  but  this  value  does  not  represent 
a node,  since  it  is  not  an  integer.  With  a similar  proof  as  the  previous  case, 
we  can  show  that  the  optimal  root  would  be  one  of  the  integers  closest  to  r„j. 
The  minimal  cost  would  be  min(cj^r„,j , cp,,,] ).  It  is  not  necessary,  however,  to 
compute  these  two  values  to  find  the  optimal  solution.  We  can  use  the  values 
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of  ri, . . . , Tp.  Let  n/,  ne  and  nh  be  the  number  of  multicast  groups  that  have 
midpoints  less  than,  equal  to,  or  greater  than  r^,  respectively.  If  [r^l  is  chosen 
as  the  root,  then  nl  groups  will  increase  by  0.5  units  each  the  total  cost,  while 
nh  groups  will  decrease  it  by  0.5  each.  The  inverse  holds  if  [r^J  is  selected. 
The  groups  with  midpoints  equal  to  will  increase  the  cost  by  0.5  anyway. 
Thus, 

C\r^\  = .5nh  - ,5nl  + .5ne  + Cr^ 


and 


c\r,r^]  - -5n/  - ,5nh  + .5ne  + 


Subtracting  these  two  equations,  we  get  that 


C[r„j  - c^r^]  =nh-nl 


and  therefore, 


nh  > nl 


Thus  the  algorithm  computes  the  optimal  solution. 

• Case  3:  p is  even,  and  there  exists  at  least  one  integer  x such  that  < x < 
rm+\.  It  can  be  shown  that  every  integer  between  and  is  an  optimal 

solution  when  p is  even.  We  will  show  that  in  particular,  is  optimal.  The 
proof  is  similar  to  the  proof  of  case  1. 
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Consider  first  the  cost  of  a node  y such  that  y < [r^].  Note  that  y < r^,  and 
thus,  as  in  case  1,  the  number  of  groups  with  mean  r,  ^ is  more  than  50%, 
or  p/2.  For  this  groups  the  cost  of  choosing  y will  decrease  by  at  most  — 2/1, 
while  the  cost  for  more  than  50%  of  the  groups  will  decrease  by  exactly  this 
amount,  increasing  the  overall  cost. 

For  y > fr„i],  the  number  of  groups  with  midpoint  r,  < y is  at  least  50%  (may 
be  more,  if  r„+i  = r„,).  Thus,  the  cost  of  selecting  y will  increase  by  at  least 
half  the  number  of  groups,  and  never  can  decrease  by  that  much,  so  Cy  could 
never  be  lower  than  Cr^.  Therefore,  [r„,]  is  optimal. 

• Case  4'  p is  even,  r„,  is  not  an  integer,  and  there  is  no  integer  between  r„  and 
This  implies  that  This  case  is  essentially  the  same  as  the  one 

discussed  in  case  1.  □ 

The  algorithm  computes  the  optimal  root  with  time  complexity  0{plogp).  Essen- 
tially, the  algorithm  computes  the  median  of  the  group  midpoints. 

Examvle  i.2  : 

Let  = 8,  p = 3,  and  the  multicast  groups  be  0,4,5,  and  2,3,4,  and  5,7.  Then 
ri  = 2.5,  T2  = 3,  and  = 6.  The  optimal  root  is  3,  and  C3  = 3 + l-f-4  = 8.  The  cost 


of  the  nodes  here  is: 
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node 

cost 

0 

16=5-1-4-1-7 

1 

13=4-l-3-h6 

2 

10=3+2-1-5 

3 

8=3+l+4 

4 

9=4+2+3 

5 

10=5+3+2 

6 

11=6+4+1 

7 

14=7+5+2 

8 

17=8+6+3 

4.3  Primary  Destinations 

The  following  algorithm  computes  the  primary  destinations  of  p multiple  groups 
for  the  linear  case  in  0{p).  The  root  of  the  multicast  tree  is  R.  The  primary  desti- 
nation for  group  i is  stored  in  PDi. 


Algorithm  i.2  Primary  destinations  for  linear  multicast  trees 
begin 

for  i=l  to  p do 

if  (/.  < R)  and  (/i,  > R)  then  PDi  = R 
else  if  /,  < R then  PDi  = hi 
else  PDi  = li 


end; 
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4.4  Minimizing  the  Real  Cost  of  the  OMT 

In  the  previous  sections,  we  have  minimized  the  sum  of  the  distances  from  the 
root  to  the  node  in  each  group  farthest  from  the  root.  The  real  cost  of  a multicast 
tree  is,  however,  the  sum  of  the  distances  from  the  primary  destination  of  each  group 
to  the  node  farther  from  it  in  that  group.  The  following  algorithm,  which  is  0{Np), 
finds  the  root  for  the  optimal  multicast  tree  in  this  sense: 

Algorithm  i.3  Optimal  root  in  terms  of  the  primary  destinations 
begin 

for  r=l  to  N-2  do  = 0; 
for  r=l  to  N-2  do 
for  i=l  to  p do 

if  {hi  < r)  or  (/,  > r)  then  Cr  = tv  + /i,  — /, 
else  Cr  = Cr  + max{r  — hi  — r}; 
r — ^^^l<r<yV— 2^r 

end; 

The  primary  destinations  can  be  computed  as  before  in  0{p)  time. 

4.5  Multiple  Trees 

As  mentioned  before,  the  multiple  group  ordering  property  can  be  satisfied  with- 
out constructing  a single  global  tree,  as  long  as  overlapping  nodes  belonging  to  two 
or  more  groups  belong  to  common  subtrees  in  the  trees  for  each  of  these  multictist 
groups.  Given  a series  of  multicast  groups,  the  graph  ^(V',  E)  is  defined  as  follows:  V 
is  the  set  of  nodes  in  the  multicast  groups,  and  for  any  two  nodes  u,  v,  (u,  v)  ^ E iff 
{u,u}  C gi  0 gj  for  any  two  groups  and  gj.  That  is,  two  nodes  are  connected  in  G 
if  and  only  if  they  both  belong  to  the  intersection  of  at  least  two  groups.  Two  nodes 
connected  by  an  edge  in  G belong  to  at  least  two  groups  i and  j.  These  two  nodes 
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must  receive  messages  addressed  to  groups  i and  j in  the  same  relative  order.  Thus, 
a subtree  containing  both  of  these  nodes  must  be  common  to  the  multicast  trees  built 
for  the  two  groups.  But  two  nodes  need  not  be  directly  connected  by  an  edge  to  re- 
quire that  they  belong  to  a common  subtree  in  the  trees  for  two  or  more  groups.  For 
example,  consider  the  following  groups:  A = {a,b,c},B  - {a,b},C  = {6, c}.  In  this 
case,  the  graph  G does  not  include  a direct  link  between  nodes  a and  c.  However,  to 
guarantee  the  multiple  group  ordering  property  it  is  required  to  construct  a subtree 
containing  all  the  three  nodes.  To  guarantee  the  ordering  property,  subtrees  for  each 
connected  component  of  G are  built.  Each  of  these  subtrees,  in  part  or  completely, 
will  belong  to  the  multicast  tree  of  two  or  more  groups.  The  subtrees  are  built  in 
the  same  way  as  we  built  multicast  trees  for  single  groups.  Once  the  subtrees  have 
been  built,  the  three  steps  described  in  chapter  3 to  construct  the  multiccist  tree  T, 
for  group  i are  executed.  These  steps  were: 

1.  First  we  need  to  determine  which  subtrees  Sj  need  to  be  part  of  T,.  Every  leaf 
in  Sj  must  belong  to  group  i.  The  nodes  and  links  in  the  subtree  that  do  not 
have  descendants  in  group  i are  not  included  in  the  subtree  for  T,.  Let  A be  the 
set  of  nodes  to  be  connected  in  step  3.  Add  the  root  of  Sj  to  set  A,  and  assign 
a depth  to  this  node  equal  to  the  depth  of  Sj.  In  case  this  node  was  assigned 
a previous  depth  x > 0,  then  assign  to  it  max{x,  depth  of  S'j}. 

2.  Consider  the  nodes  j in  group  i that  did  not  belong  to  any  connected  compo- 
nent. If  node  ] can  be  covered  by  any  of  the  subtrees  found  in  the  previous 
step  without  increasing  the  cost  of  the  subtree,  then  considered  this  node  as 
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covered  by  the  subtree  (this  may  affect  the  traffic  generated  by  the  subtree,  but 
not  its  cost).  This  node  is  inserted  below  the  closest  node  to  it  already  in  the 
subtree  that  keeps  the  subtree  a shortest-path  tree.  The  longest  possible  path 
already  in  the  tree  is  used  as  a shortest-path  from  the  root  to  node  j.  This 
helps  in  reducing  the  traffic.  The  other  nodes  that  cannot  be  inserted  without 
increcising  the  cost  of  the  subtree  are  added  to  A and  a depth  of  0 is  cissigned 
to  them. 

3.  Connect  the  nodes  in  set  A. 

We  discuss  now  how  the  nodes  in  A are  connected  for  linear  networks.  Note 
that  now  the  cost  of  node  i is  given  by  c,  = ma.x d{i,j)  -f  depth(j).  First  we 
find  the  smallest  node  a E A such  that  depth{a)  > a — x + depth[x)  for  any  node 
X ^ A with  X < a.  Similarly,  node  b E A is  the  largest  node  in  A such  that 
depth{b)  > y -b-\-  depth{y)  for  any  node  y E A with  y > b.  Note  that  it  may  happen 
that  a = b.  Let  r = (a  -|-  6 + depth{b)  — depth{a))f2.  If  r'  < a then  r = a.  Else  if 
r > b then  r — b.  Else  r = r' . 

Lemma  J.l  r is  the  optimal  root. 

Proof: 

Assume,  for  simplicity,  that  r'  is  an  integer.  If  r'  < a or  r'  > fe  then  r is  clearly  the 
optimal  root.  Otherwise  r lies  exactly  in  the  middle  of  a and  b + depth{b)  - depth{a), 
and  therefore  r — a -f  depth(a)  = b — r + depth[b)  = Cr.  Let  x be  an  optimal  root.  If 
X < r then  > b — r + depth{b)  = c^,  which  is  a contradiction.  Similarly,  if  x > r, 
then  Cx  > Cr.  Thus  x = r,  and  hence  the  proof  is  complete.  □ 
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r 


Figure  4.1.  Form  of  each  tree  in  a linear  network 


The  final  tree  T,  for  group  i must  have  this  form  given  in  figure  4.1,  where  r is  the 
optimal  root,  and  / and  h are  the  lowest  and  highest  nodes  in  group  z,  respectively. 

Some  of  the  nodes  in  the  path  from  r to  I may  have  a linear  portion  to  their  right 
side.  Similarly,  some  of  the  nodes  in  the  path  from  r to  h may  have  a linear  portion 
to  their  left  side.  These  nodes  were  roots  of  the  subtrees  built  for  the  connected 
components  of  graph  G.  The  nodes  in  the  other  side  of  these  subtrees  can  be  reached 
by  the  paths  from  r to  I and  h,  preserving  the  structure  of  the  subtrees.  This  helps 
reduce  the  traffic  generated  by  the  multicasts. 


Example  4-3  : 


Let  the  multicast  groups  be  given  by:  A = {0,2,3},.S  = {2, 3, 4, 5, 7},  C = {3,5},Z)  = 
{4, 6,7,8}. 

The  connected  components  are  given  by:  {2,3,5}  and  {4,7}.  The  subtrees  for 

these  two  connected  components  are,  respectively: 

3 5 


2 


4 


4 6 


5 
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For  group  A,  only  the  first  subtree  must  be  part  of  its  multicast  tree.  We  need  to 
connect  node  0 (depth=0)  and  node  3 (depth=l).  The  optimal  root  is  node  2,  and 
the  multicast  tree  is  given  by: 


2 


1 3 

0 2 


For  group  B we  need  to  connect  both  subtrees,  having  both  roots  3 and  5 a cost 
of  2.  The  tree  for  group  B is: 


For  group  C only  a segment  of  the  first  subtree  is  needed.  For  group  D,  only  the 
second  subtree  is  needed.  Also,  when  nodes  6 and  8 are  considered  in  step  3,  node  6 
is  discarded  since  it  is  covered  by  the  subtree,  and  node  8 need  to  be  connected  with 
node  5.  The  trees  for  groups  C and  D are  shown  below: 


3 

4 


6 


8 


5 


7 
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The  primary  destinations  for  groups  A,  B,  C,  D are  respectively,  nodes  2, 4, 3,  and 
6.  In  this  case  the  node  that  is  primary  destination  a maximum  number  of  times  is 
primary  destination  for  25%  of  the  multicast  groups.  □ 

In  the  previous  example,  the  cost  to  complete  the  multicast  is  given  by  2+3+2-f-3=10, 
while  the  traffic  generated  is  4+7+2+6=19.  Using  the  single  global  tree  approach 
for  this  example  we  would  have  obtained  the  multicast  given  by: 


4 


3 5 


2 6 

1 7 

0 8 

with  a cost  of  3+3+1+4=11,  and  a traffic  of  3+5+2+4=14.  Node  4 is  now  primary 
destination  for  75%  of  the  groups. 

The  following  graphs  illustrate  our  results  from  simulations  ran  for  a linear  net- 
work with  64  nodes.  To  increase  the  number  of  nodes  in  the  intersections  of  the 
multicast  groups,  we  assumed  the  existence  of  ’’hot  nodes”  in  the  network.  In  our 
simulations,  the  first  25%  of  the  nodes  in  the  network  have  a fixed  probability  of  0.2 
of  being  included  in  a multicast  group,  while  the  rest  have  only  probability  of  0.1. 
Figure  4.2  shows  the  cost  of  the  multicasts  (assuming  that  one  multicast  message  is 
addressed  to  each  multicast  group)  for  the  randomly  generated  multicast  groups.  The 
graph  shows  that  the  cost  of  the  multicasts  is  lower  for  the  multiple  trees  model  than 
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for  the  single  global  tree  model.  Figure  4.3  shows  that  the  traffic  generated  is  higher 
for  the  multiple  trees  approach.  As  mentioned  before,  this  result  Wcis  expected,  since 
in  this  model  the  same  message  may  be  required  to  be  sent  over  the  same  link  (in 
the  same  or  in  opposite  directions)  more  than  once.  The  last  example  given  above 
illustrates  this  fact.  Thus  the  traffic  was  expected  to  be  higher  in  the  multiple  trees 
approach  to  guarantee  the  multiple  group  ordering  property.  The  third  graph  shows 
what  is  probably  the  most  significant  advantage  of  the  multiple  trees  model:  much 
higher  variability  in  the  primary  destinations.  In  the  single  global  tree  model,  the 
same  node  (the  root  of  the  tree)  must  be  the  primary  destination  for  most  of  the 
multicast  groups.  In  that  case  we  do  not  deviate  too  much  from  a centralized  solu- 
tion to  the  multiple  group  ordering  problem,  which  has  the  potential  disadvantage  of 
developing  a bottleneck  in  the  root.  In  the  multiple  trees  approach  the  function  of 
being  a primary  destination  is  much  more  distributed  among  different  nodes  in  the 


network. 


Cost  of  Multicast 
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Number  of  Groups 


Figure  4.2.  Cost  for  Linear  Network  (N=64) 


Traffic 
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Number  of  Groups 


Figure  4.3.  Traffic  for  Linear  Network  (N=64) 


% of  Identical  Primary  Dest 
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Number  of  Groups 


Figure  4.4.  Primary  Destinations  for  Linear  Netwrok  (N=64) 


CHAPTER  5 

OMT  FOR  RING  NETWORKS 
5.1  Single  Group 

Let  Rn  be  the  ring  network  with  nodes  0,1,2,...,N-1.  Let  D be  the  destination 
list.  A path  from  node  x to  node  y is  denoted  by  (x,y).  Nodes  x and  y are  the 
endpoints  of  this  path.  All  other  nodes  in  the  path  are  called  interior  points.  The 
set  of  interior  points  for  a path  P is  denoted  by  I{P).  For  every  pair  of  nodes  x 
and  y there  are  two  paths  in  that  connects  them.  When  it  becomes  necessary  to 
differentiate  between  them,  extra  nodes  in  the  path  would  be  added  to  denote  it,  like 
for  example,  (x,x  + 1, ...,  y)  is  the  path  from  x to  y that  goes  through  node  x + 1.  For 
simplicity,  sometimes  the  set  D is  represented  without  brackets  and  without  commeis 
between  its  elements.  For  example,  D = 024  denotes  the  multiccist  group  composed 
of  the  nodes  0,  2 and  4.  The  complement  of  a path  P (x,  x + 1, ...,  y)  is  defined  a5  the 
path  (x,x  - l,...,y),  and  viceversa.  It  is  denoted  by  P.  Here,  additions  are  modulo 
N. 

Let  P be  the  longest  path  in  the  ring  such  that  if  x € I{P)  then  x ^ D. 

Lemma  5.1  The  optimal  root  is  the  midpoint  of  the  complement  of  path  P. 

Proof: 

Let  m be  the  midpoint  of  P.  Let  a and  b be  the  endpoints  of  P.  Clearly,  a and  b are 
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elements  of  D.  Assume  that  if  the  number  of  nodes  in  path  P is  even,  then  we  define 
m to  be  the  node  that  is  closer  to  b.  Thus,  = d{m,a)  < d(m,b)  + 1.  We  need 
to  show  that  m is  an  optimal  solution,  that  is,  for  any  other  node  x in  the  network, 
Cm  < Cx-  It  should  be  clear  that  any  node  in  P satisfies  this  condition.  We  need  to 
show  that  it  is  also  true  for  nodes  in  I{P).  Let  x € I{P)- 

Case  1:  m ^ D.  Clearly,  c„,  = d{m,a)  < d{m,x).  Therefore,  Cm  < c^,  and  we  are 
done.  Note  that  a is  the  node  in  D farthest  from  m. 

Case  2:  m ^ D.  Assume,  by  contradiction,  that  < Cm-  Let  y and  z be  the 
elements  in  D closest  to  m such  that  y is  closer  to  a and  z is  closer  to  b.  Note  that 
it  may  be  the  Ccise  that  a = y or  b = z.  See  the  figure  below  to  locate  the  nodes  we 
have  defined  so  far. 


P 


P 


We  have  that 


d{m,a)  = d{m,y)  + d{a,y) 


d{x,y)  = d{x,a)  + d{a,y) 
Subtracting  the  above  equations,  we  get  that 


d{m,a)  - d{x,y)  = d{m,y)  - d{x,a) 
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Since  Cl  < Cm,  we  have  that  d(m,y)  — d{x,a)  > 0.  This  implies  that  d(m,y)  > 
<f(x,a).  Similarly,  we  can  show  that  d{m,z)  > d(x,b).  Combining  these  two  inequal- 
ities, we  get  that  the  path  from  y to  z that  goes  through  m is  greater  than  the  path 
from  a to  b that  goes  through  x.  But  this  contradicts  the  maximality  of  P,  since  by 
construction  the  first  path  does  not  include  any  nodes  in  D.  Thus  the  lemma  has  to 
be  true.  □ 

The  complement  of  path  P,  and  thus  the  optimal  root,  can  be  found  with  an  0(N) 
algorithm. 

Example  5.1  : 

Let  N = 8.  Let  D = 1,2, 5 (shown  in  dark  in  the  figure). 


0 


4 


P is  (1,0, 7, 6,5),  and  P = (1,2, 3, 4, 5).  The  midpoint  of  P is  node  3,  and  this  is 
the  optimal  root  with  C3  = 2.  □ 

5.2  Multiple  Groups 

Let  n be  the  maximum  distance  possible  from  any  root  to  any  node  in  a shortest- 
path  tree.  In  a ring,  n = \_Nj2\ . Any  multicast  tree  must  have  the  same  shape  (see 
figure  below).  For  a given  root  r,  and  a shortest-path  tree  rooted  at  r,  let  L,>  be  the 
maximum  distance  from  r to  the  last  node  of  group  i in  the  left  side  of  the  tree,  and 
Rtr  the  corresponding  distance  in  the  right  side,  as  shown  in  figure  5.1. 
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Figure  5.1.  Form  of  tree  in  a ring  network 

Clearly,  0 < Lir,  Rir  < n for  all  i.  Before  presenting  the  algorithm,  let  us  see  first 
how  and  why  it  works,  as  we  describe  the  purpose  of  the  variables  used  in  it. 

The  algorithm  first  computes  (going  through  the  elements  in  D)  the  values  of  Rio 
and  Lio  for  all  the  groups.  In  general,  once  Rij  and  Z/,j  are  computed,  to  compute 
the  values  of  Rij+\  and  Lij^i,  since  we  are  moving  clockwise  through  the  ring,  we 
may  need  to  subtract  1 from  Rij  and  add  1 to  L{j.  The  node  farthest  from  node  j 
to  the  right  is  now  one  unit  closer  to  node  j + 1,  and  the  node  farthest  from  j to 
the  left  is  now  one  unit  farther  from  j T 1.  There  are  some  special  cases,  however, 
that  complicate  things  a little  bit.  The  highest  value  that  can  get  is  n,  so  when 
Lij  = n vfe  cannot  increase  it  to  compute  T,j+i.  Also,  the  lowest  value  can  get  is 
0,  so  again  we  need  to  compute  Rij+i  in  some  other  way  when  Rij  = 0.  It  is  necessary 
to  deal  with  these  cases  separately.  There  is  one  more  special  case.  This  one  occurs 
when  L,j  = 0. 

When  Lij  — 0,  and  we  move  to  node  j + 1,  Lij+i  does  not  necessarily  gets  increased, 
but  instead  keeps  its  value  of  0.  Consider,  for  example,  the  ring  with  5 nodes.  If 
D={1,2},  Lio  = 0,  Rio  = 2,  but  Lu  remains  at  0.  while  Ru  decreases  to  1.  This 
situation  occurs  whenever  all  the  nodes  in  D are  covered  by  the  right  portion  of  the 
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tree,  and  j ^ D.  The  way  to  solve  this  case  is  to  keep  a counter  that  stores  the 
distance  between  node  j and  the  closest  element  to  j in  D.  Essentially,  the  counter 
keeps  track  of  how  many  iterations  (when  L = 0)  we  would  have  to  keep  L = 0.  This 
counter  is  called  counter _L0[i]  in  the  algorithm,  and  it  gets  initialized  every  time  L;j 
becomes  0,  with  the  distance  between  node  j + 1 and  the  node  following  j + 1 in  D. 

When  Rij  = 0,  it  means  that  the  destination  list  is  covered  completely  by  the  left 
side  of  the  multicast  tree  when  j is  the  root.  It  will  remain  at  0 until  L becomes  n. 

Now,  assume  that  Lij  = n.  When  j + 1 is  the  root,  then  the  node  farthest  from 
node  j counterclockwise  (let  us  call  this  node  j)  would  be  covered  by  the  right  side 
of  the  tree.  If  N is  odd,  the  distance  from  j + 1 to  j is  n,  while  when  N is  even  it 
is  n — 1.  This  value  is  assigned  to  Rij+i-  Now,  Z/,_,+i  would  be  the  distance  between 
node  j + l and  the  node  following  j in  Z),  ji  + 1,  if  this  node  is  not  covered  by  the  right 
side  of  the  tree,  which  is  the  path  {j  + !,_;  + 2, . . . ,j).  Otherwise,  the  destination  list 
is  covered  completely  by  this  path,  and  L,j+i  becomes  0.  The  value  of  J is  computed 
by  j = {N  + j — n)  mod  N.  Let  CWi{a)  be  the  node  in  group  i closest  to  node  a 
clockwise,  including  a itself. 

Algorithm  5.1  Ring-Multiple  Groups 
begin 

1.  for  i=l  to  p do 

compute  /?,o,  T,o; 

if  L,o  = 0 then  counter  J^0[i]  = ClT,(0) 
else  counter _L0[i]  = 0; 

2.  for  i=l  to  p do 

for  i=0  to  N-2  do 
if  Lij  = n then 

if  N is  even  then  R,j+i  = n — 1 
else  Rij+i  = n; 


z = CWiU  + 1); 
if  X ^ path  {j  + l,j -1-2,...,;) 
then  Lij+i  = d{j  -|-  l,x) 
else 

Lij+l  = 0, 

counter dj0[i]  = d{j  -f  1,CW,(;  4-  1)) 

else 

if  counter_L0[i]  > 0 then 

Rij^i  = lij  Ij 
•f'tj+i  — 0, 

counter  _L0[i]  = counter  .L0[i]  — 1; 

else 

if  Rij  / 0 then  Rij+i  = Rij  — 1 else  = 0 

Rij+l  — Rij  4"  1) 

3.  select  an  optimal  root  j that  minimizes 
Ef=i  max  {Lij,  Rij] 

end; 


Example  5.2  : 

Consider  the  following  ring,  with  D = 2,3, 7. 


The  algorithm  computes  the  following  values  for  R and  L: 
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Node  j 

Rj 

0 

3 

6 

1 

6 

0 

2 

5 

0 

3 

4 

1 

4 

3 

2 

5 

2 

3 

6 

1 

4 

7 

0 

5 

8 

0 

6 

9 

6 

6 

10 

6 

3 

11 

5 

4 

12 

4 

5 

Nodes  4 and  5 are  both  optimal  with  cost  3.  □ 

The  algorithm  requires  then  0{Np)  time  complexity.  This  is  better  than  the 
0{Nm)  required  to  compute  the  cost  of  every  node  going  through  all  the  elements 
of  the  destination  lists  for  every  node.  Note  that  m itself  may  be  in  the  order  of  Np. 
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5.3  Primary  Destinations 

The  primary  destinations  for  the  multicast  tree  constructed  for  rings  can  be  com- 
puted again  in  0{p).  Let  and  Rij  be  the  values  computed  by  Algorithm  5.1  for 
group  i and  node  j.  The  optimal  root  found  by  the  algorithm  is  the  node  j that 
minimizes  max{L,j,  For  each  group  i,  the  primary  destination  PD{  is 

computed  as  follows: 

if  (L,j  > 0)  and  {Rij  > 0)  then  PDi  — root 
else  if  Lij  = 0 then  PDi  = CWi{root) 
else  PDi  = CCWi{root) 

where  CCWi{root)  is  the  closest  node  to  the  root  counterclockwise. 

5.4  Minimizing  the  Real  Cost  of  the  OMT 

For  any  shortest-path  tree  the  sum  of  the  distances  from  the  primary  destinations 
to  the  nodes  farther  from  them  in  each  group  can  be  minimized  in  the  following  way. 
Instead  of  finding  the  value  that  minimizes  max{L,j,  we  find  the  node  r 
that  minimizes  the  following  expression: 


E 


«=i 


max{Lij^Rij]  if  Lij,  Rij  > 0 

Rij  — d{r,x)  if  Lij  = 0 (x  = CkF,(r)) 

Lij  - d{r,  x)  if  R,j  = 0 (x  = CCWi{r)) 


From  all  the  shortest  path  trees,  this  method  will  give  us  the  optimal  root  in 
terms  of  the  distances  from  the  primary  destinations.  However,  using  this  method 
we  always  produce  a shortest  path  tree,  and  the  optimal  multicast  tree  does  not 
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necessarily  have  this  property.  Consider,  for  example,  the  ring  with  N = 7,  and 

groups  016,  23,  34,  and  45.  The  OMT  in  this  case  is  given  by: 

0 0 

x\ 

} 6 


5 

4 

3 


/\ 
1 ( 


None  of  these  trees  is  a shortest  path  tree,  and  thus  the  method  given  above  will 
not  construct  the  OMT.  Let  Rijk  and  be  the  cost  of  the  right  and  left  side  of  the 
multicast  tree  for  group  i when  node  j is  the  root,  and  exactly  k elements  of  group 
i are  covered  by  the  right  side  of  the  tree.  Note  that  now  the  tree  does  not  have  to 
be  a shortest-path  tree.  Let  y)  and  be  the  distances  between  nodes 

X and  y counterclockwise  and  clockwise,  respectively.  Also,  let  CWi{a)  (respectively 
CCWi{a))  be  the  node  in  group  i closest  to  node  a clockwise,  including  node  a.  The 
following  algorithm,  which  is  0{Nmp),  constructs  the  OMT  for  any  ring  network, 
and  finds  the  primary  destinations.  Remember  that  m = m,,  and  m,  is  the 

number  of  nodes  in  group  i. 

The  values  of  R{jk  and  Lijk  are  computed  in  the  following  way.  When  k = 0, 
all  nodes  in  group  i must  be  covered  by  the  left  side  of  the  tree,  so  R.j*.  = 0,  and 
Lijk  is  the  counterclockwise  distance  between  node  j and  the  closest  element  in  the 
group  to  node  j + 1 clockwise;  for  other  values  of  k^  Rijk  is  the  clockwise  distance 
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between  j and  x,  where  x is  the  closest  node  clockwise  to  the  one  used  in  the  previous 
iteration  when  k was  ^ — 1,  while  L{jk  has  to  cover  all  nodes  following  node  x in  the 
group.  When  k = m,,  then  all  nodes  are  covered  by  the  right  side  of  the  tree,  and 
thus  Lijk  = 0.  If  node  j is  an  element  of  group  i,  then  only  m,  — 1 nodes  have  to  be 
covered  by  the  tree,  and  x will  be  j again  when  /:  = m,  — 1.  In  this  case  we  set  Lij^, 
and  Rijmi  to  their  previous  values  and  k is  increased  to  exit  the  loop. 


Algorithm  5.2  OMT  for  ring  networks 


begin 

for  i=l  to  p do 
for  j=0  to  N-1  do 


x=j; 

for  k=0  to  m,  do 
R,jk  = d^{j,  i); 
i(  k = TUi  then  Lijk  = 0 
else  L,,k  = d^^{j,CW,{x+l))- 
X = CW,{x  + 1); 
if  X = J then 
Rijk+l  — djijk] 

R-ijk+1  — Rfjki 

k=k+l; 


COSTi.k  = 


end; 


mox^Lijk^  Rijk} 
R.jk-d{j,CW,{j)) 

L„k-d{j,CCW,{j)) 


if  Lijk,  Rijk  > 0 
if  Lijk  - 0 
if  R,jk  = 0 


The  OMT  can  be  found  minimizing  COST^k  = Ef=i  COSTijk-  The  optimal  root 
and  the  routing  are  given  by  the  indexes  j and  k that  minimize  COSTjk- 


5.5  Multiple  Trees 

The  three  steps  as  discussed  in  section  4.5  to  build  the  multicast  tree  T,  for  group 
i once  the  subtrees  for  the  connected  components  of  the  graph  G as  defined  have 
been  built  can  be  used,  as  mentioned  before,  by  the  different  network  topologies  we 
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consider.  We  need  to  discuss  how  to  perform  the  third  step  of  the  algorithm,  that  is, 
how  to  connect  the  nodes  in  set  A. 

First  we  compute  the  optimal  root  r'  for  the  set  A assuming  that  all  nodes  in 
A have  depth  0.  To  do  this  the  algorithm  described  in  section  5.1  to  compute  the 
optimal  root  for  a single  group  can  be  used.  Now  we  determine  the  nodes  a,b  £ A 
with  the  following  properties: 

• distance{r  ,a)  = clockwise  distance  (r’,a) 

• distance{r\a)  + depth{a)  is  maximum  for  all  the  nodes  in  A satisfying  the 
previous  property 

• distance{r  ,b)  = counterclockwise  distance  {r',b) 

• distance{r  ,b)  + depth{b)  is  maximum  for  all  the  nodes  in  A satisfying  the 
previous  property 

Taking  r as  the  root,  let  x and  y be  the  costs  of  the  left  (clockwise)  and  right 
(counterclockwise)  sides  of  the  multicast  tree.  The  cost  of  the  root  node  r is  thus 
ma.x{distance{r' ,a)  + cost{a),distance{r\b)  + depth{b)]  = max{x,y}.  Essentially, 
to  compute  the  optimal  root,  we  need  to  move  \ x - y \ j2  nodes  in  the  direction  of 
node  a if  I > y or  in  the  direction  of  6 if  y > x.  Then  we  would  obtain  the  midpoint 
r between  distance{r\a)  -f  depth{a)  and  distance{r' ,b)  + depth(b).  The  cost  of  the 
left  side  of  the  tree  would  be  as  close  as  possible  to  the  cost  of  the  right  side.  Clearly, 
this  would  yield  the  optimal  root. 
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The  final  tree  has  the  same  form  cis  the  tree  for  linear  networks  shown  in  figure 
4.1. 

Simulations  ran  for  the  ring  network  with  64  nodes  under  the  same  conditions 
ais  the  ones  ran  for  linear  networks  show  essentially  the  same  results.  The  cost  of 
the  multicasts  is  better  for  the  multiple  trees  approach,  while  at  the  same  time  it 
generates  more  traffic.  The  primary  destinations,  as  for  linear  networks,  have  much 
more  variability  under  the  multiple  trees  approach. 


Cost  of  Multicast 
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Number  of  Groups 


Figure  5.2.  Cost  for  Ring  Network  (N=64) 


Traffic 
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Figure  5.3.  Traffic  for  Ring  Netwrok  (N=64) 


% of  Identical  Primary  Dest. 
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Number  of  Groups 


Figure  5.4.  Primary  Destinations  for  Ring  Netwrok  (N=64) 


CHAPTER  6 

MULTICAST  TREES  FOR  MESH  NETWORKS 
6.1  Single  Group 

Let  M be  an  m*n  mesh.  Let  D be  a set  of  points  from  the  grid.  The  cost  of  node 
i is,  as  before,  defined  by  c,  = maxjgo  d(f,  7).  The  distance  in  a grid  between  the 
points  (xi,X2)  and  (j/i,j/2)  is  given  by  | xj  — ?/i  | + | X2  — 1/2  |-  The  first  thing  we 
need  to  do  is  find  the  node  r with  minimum  cost. 

Lemma  6.1  Let  y be  an  optimal  solution.  For  any  node  x,  if  Cy  < Cx,  then  x has  a 
neighbor  or  a diagonal  z such  that  < c^. 

Proof: 

Let  X = (1,7),  y = {h,l).  Assume  that  Cy  < c^.  This  implies  that  x ^ y.  Assume, 
w.l.g.,  that  i < k.  Let  max„g£)  d(a,  n)  = d{a,ma),  For  any  node  a,  rria  is  the  node  in 
D farthest  from  node  a.'  Since  y is  optimal,  we  know  that  d{y,my)  < d(x,mj;),  and 
d{y,my)  < d{z,m^)  for  any  neighbor  z. 

Case  1:  x and  y have  the  same  second  component  (7  = /).  Then  choose  z = 
(i  + 1,7).  We  need  to  show  that  < Cx.  The  costs  of  z and  x are  given  by 
Cz  = d{z,m^)  and  Cx  = d{x,mx).  We  know  that  y is  at  leeist  as  close  to  as  z is 


‘Note  that  rria  may  not  be  unique. 
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j 


i 

i+1 

k 


Figure  6.1. 

(note  that  maybe  y = z).  Thus,  the  first  component  of  is  at  least  i + 1.  Given  x 
and  y as  shown  in  figure  6.1,  must  be  somewhere  in  the  shaded  area. 

Since  X and  z have  the  same  second  component,  it  is  clear  that  d(z,m;i)  < 
d{x,m^)  < d{x,mx).  Thus,  c*  < c^. 

Case  2:  x and  y have  different  second  component.  Assume,  w.l.g.,  that  / > j. 
That  is,  y is  located  below  and  to  the  right  of  x in  the  grid.  Let  z\  = (i  + l,j),Z2  = 
{ij  -f  1),23  = {i  + l,j  + 1).  We  need  to  show  that  c^,  < for  i=l,2,  or  3.  We  prove 
this  result  by  contradiction. 

Assume  that  d{x,  m^)  < d(z,,  m^, ) Vi  = 1 , 2, 3.  We  know  that  d(x,  m^, ) < d{x,  rrix)  < 
d(2i,mj, ).  That  is,  x is  closer  to  than  Z\  is.  This  implies  that  m^,  has  first  com- 
ponent < i.  A similar  relation  for  shows  that  it  has  second  component  < j. 
The  fact  that  d^y^m^^)  < d(^Zi,m2,)  for  i = l,2  implies  that  m^,  has  second  component 
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i 

i+1 

k 


> j + 1,  and  that  has  first  component  > i + 1.  Also,  d{x,mx)  < d{z3,m^^), 
which  implies  that  nix  has  first  component  < i or  second  component  < j,  but  not 
both  (since  d{y,nix)  < d{x,nix)).  Given  x and  y as  shown  in  figure  6.2,  and 
must  be  in  the  shaded  areas,  while  can  be  in  either  of  them. 

Now,  consider  a shortest  path  from  m^,  to  x,  and  from  x to  This  path  has 
to  be  a shortest  path  from  m^,  to  m^j.  x must  be  the  midpoint  of  this  path.  That 
is,  = d(m^,,x)  + d(x,m^J,  and  | d(x,m^, ) -d(x,m^J  |<  1.  This  last 

inequality  must  hold,  or  otherwise,  if  for  example  d(x,  ) > d(x,  ) + 1,  then  m^, 
would  be  the  point  farther  from  Z2,  and  not  as  it  is. 


j j+i 


X 

m 

"3 

y 

Figure  6.2. 
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Since  X and  2,  are  adjacent,  and  d{x,mx)  < d{zi,m2-)  for  i=l,2,  this  implies 
that  d{x,mx)  < d{zi,m.z.)  for  i=l,2.  ^ Thus  it  must  be  the  Ccise  that  d{x,mx)  = 
ma,x(d{x,m^^),d{x,m^^)).  Assume,  w.l.g.,  that  d{x,m^^)  > d{x,m^^).  d{y,mx)  < 
d{x,mx)  implies  that  d(i/,m^, ) < d(i,mz, ).  Therefore,  d(y,rrix^)  > d(x,m^^)  > 
d{x,mx)  — 1.^  Thus,  d{y,m.2^)  > d{x,mx).  This  implies  that  Cy  > Cx,  which  is  a 
contradiction.  Therefore,  the  statement  of  the  lemma  must  be  true.  □ 

The  lemma  states  that  if  a node  is  not  an  optimal  root,  then  it  must  have  a 
neighbor  or  a diagonal  with  lower  cost.  The  idea  then  is  to  start  with  a good  candidate 
for  an  optimal  solution,  and  then  if  it  turns  to  be  not  an  optimal  solution,  we  move 
in  the  direction  of  an  optimal  until  we  eventually  converge  to  one.  We  select  as  the 
initial  choice  the  node  such  that  the  first  and  second  components  are  the  midpoints  of 
the  lists  of  first  and  second  components,  respectively,  of  the  elements  in  the  multicast 
group.  Then  we  check  if  that  node  has  a neighbor  with  less  cost.  If  it  doesn’t,  then 
it  is  optimal.  Otherwise  we  move  to  that  new  node,  and  repeat,  until  we  reach  an 
optimal  solution. 

In  a linear  network,  computing  the  midpoint  of  the  smallest  node  and  the  highest 
node  gives  us  the  optimal  root.  In  a 2-dimensional  mesh,  computing  the  midpoint  that 
way  for  each  dimension  separately  to  get  the  two  coordinates  of  a point  not  necessarily 
will  produce  the  optimal  root.  The  distance  from  that  point  to  the  lowest  and  highest 
rows  and  columns  in  the  mesh  containing  elements  of  D would  be  minimized,  but 
the  sum  of  the  horizontal  and  vertical  distances  to  reach  the  farthest  node  may  not 
^Two  adjacent  nodes  cannot  have  the  same  distance  to  any  node. 

^The  first  inequality  follows  from  the  fact  that  the  path  from  to  x,  and  from  x to  is  a 
shortest  path  from  m^,  to  together  with  the  fact  that  d(y,  m,,)  < d(x,m^,). 
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be  minimized.  Let  h be  the  maximum  possible  of  nodes  to  which  we  would  have  to 
jump  in  the  algorithm  to  compute  the  optimal  root.  To  find  the  value  of  h,  a worst 
case  scenario  should  be  considered.  An  example  of  this  is  illustrated  in  the  figure 
below. 


Let  D be  as  shown  in  the  figure.  The  first  choice  for  the  root  is  the  node  (4,4). 
This  point  minimizes  the  distance  to  the  lowest  column  and  row  with  elements  from 
D.  However,  it  is  not  optimal.  Let  dj  be  the  distance  from  that  point  to  the  lowest 
(or  highest)  column  that  contains  an  element  of  D,  and  d,2  the  corresponding  distance 
to  the  lowest  (or  highest)  row.  The  maximum  cost  of  the  point  computed  is  di  + ^2, 
as  it  is  in  our  example.  From  one  side  of  the  potential  tree  the  cost  is  di  4-^2,  and 
from  the  other  side  it  is  only  di  (in  this  case).  We  would  like  to  ’’balance”  this  values. 
Moving  h nodes  to  the  right,  we  can  balance  this  distances,  making  di  h as  close 
as  possible  to  di  + ^2  — h.  This  implies  that  h = [d2/2J.  d2  is  at  most  [m/2j,  so  h 
is  0(m/4).  In  this  example,  dj  = d2  = 3,  and  h = 1.  The  algorithm  thus  has  time 
complexity  0{m  | D |)  in  the  worst  case. 


Example  6. 1 : 
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Let  M be  the  5*6  grid,  and  let  D = {(1, 3),  (2, 5),  (4, 2),  (4, 3),  (5, 1),  (5, 6)},  as 
shown  in  the  figure  below. 


1 2 3 4 5 6 


Our  first  choice  for  optimal  root  is  the  node  x = ( [(1  +5)/2j , [(1  -f-6)/2j ) = (3, 3). 
C3,3  = 5.  Now  we  check  if  any  neighbor  or  diagonal  has  a lower  cost.  In  fact,  node 
(4,3)  has  cost  4.  This  node  does  not  have  a neighbor  or  a diagonal  with  lower  cost, 
hence  it  is  optimal.  □ 

The  following  lemma  will  be  useful  for  the  solution  of  the  multiple  group  case. 
Lemma  6.2  Let  x=(i,j),  y=(k,l),  k > i,  I > j,  Cy  < c^,  and  Z3  = (t  + l,j  + 1).  Then 

Cjj  ^ Cj;  . 

Proof: 

The  only  way  that  d{x,mj;)  could  be  smaller  than  is  if  rUx  has  first  com- 

ponent < i,  and  second  component  < j.  This  cannot  be  the  case,  since  d[y,mx)  < 
d[x,mx).  Thus  the  lemma  is  true.  □ 

The  problem  now  is  to  determine  the  routing  to  construct  the  multicast  tree  using 
the  optimal  root. 

3 he  problem  of  finding  a Steiner  tree  given  a set  of  points  in  a plane  using  rectilin- 
ear distances  has  been  proved  to  be  an  NP-Complete  problem  [15].  Here  we  propose 
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a heuristic  algorithm  to  find  a suboptimal  multicast  tree  (SOMT).  We  wish  to  find 
a tree  with  the  following  properties: 

• V{T)  C V{M) 

• the  maximum  of  the  distances  from  the  root  to  every  node  in  the  destination 
list  is  as  small  as  possible 

• if  i?  is  the  root  of  the  tree,  then  the  following  holds:  dx{x,  R)  = d^/(x,  R) 
Vx  e SOMT 

• no  tree  exists  satisfying  the  three  properties  above  that  has  fewer  nodes  than 
this  one 

The  third  condition  above  implies  that  the  heuristic  algorithm  will  construct  a 
shortest-path  tree  for  the  nodes  in  the  multicast  group  D.  Constructing  a shortest- 
path  tree  with  lowest  traffic  rooted  at  node  R (grid  M in  the  figure  below)  is  equiv- 
alent to  finding  a rectilinear  directed  Steiner  tree  rooted  at  i?  in  a directed  mesh  in 

which  all  the  links  are  directed  away  from  R {M‘  in  the  figure). 

M M‘ 


R 


The  problem  of  finding  a rectilinear  directed  Steiner  tree  has  been  known  in  the 
literature  as  the  Rectilinear  Steiner  Arborescence  (RSA)  problem.  The  problem  here 
is  to  find  a Steiner  directed  tree  for  a set  D of  nodes  lying  in  the  grid,  rooted  at  point 


r' 
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R,  using  only  horizontal  and  vertical  arcs  oriented  away  from  R.  Rao  et  aJ.  [32] 
have  studied  the  Rectilinear  Steiner  Arborescence  problem  when  the  set  D lies  in  the 
first  quadrant  of  , and  the  directed  tree  is  rooted  at  the  origin.  They  proposed 
a heuristic  algorithm  in  which  a Steiner  tree  H for  D is  generated  by  iteratively 
substituting  min[p,  q)  for  pairs  of  points  p,q  'm  D until  the  origin  remains,  where 
p = (xp,yp),  q = and  min{p,q)  = {min{xp,yp),min{xq,yg)).  The  points  p 

and  q are  chosen  to  maximize  ||  min{p,q)  ||  over  all  the  points  currently  in  D.  They 
have  shown  that  1{H)  < 21{RSM A),  where  1{A)  is  the  number  of  links  in  the  tree  A 
and  RSM A is  an  RSA  A with  minimum  /(A). 

A similar  heuristic  is  presented  here  that  generates  a SOMT  for  an  mxn  grid  when 
D lies  anywhere  in  the  grid.  For  two  grid  points  p,  q,  let  maXo{p,  q)  be  the  grid  point 
that  maximizes  the  overlapping  in  the  shortest  paths  from  the  origin  to  the  points  p 
and  q.  maXo{p,  q)  is  defined  in  the  following  way,  given  that  p = (ip,  j/p),  q = (x,,  y,), 
for  z = X or  y coordinates;  if  Zp,  z,  > 0 then  the  z-coordinate  of  maxo  is  given  by 
min[zp,  Zg).  Else  if  Zp,z,  < 0 then  the  z-coordinate  is  max{zp,  Zg).  Otherwise  it  is  0. 
We  assume  that  the  root  of  the  tree  is  the  origin. 

Let  denote  the  four  lines  lying  in  the  mxn  grid  given  by  | x | -f  | y |=  z.  A set 
of  points  P is  a cover  for  a set  of  points  Q ii^q  & Q there  exists  a point  p € P with 
a path  from  p to  q.  Let  be  the  subset  of  Q lying  above  L,.  Then  let  MC{Q,z) 
denote  the  set  of  points  on  L^with  minimum  cardinality  that  covers  Q^.  It  follows 
directly  from  this  definition  that  any  RSA  for  D intersects  at  least  MC{D,z) 
times  for  every  z.  The  nodes  in  the  initial  set  D are  called  sinks. 
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A heuristic  algorithm  generates  a SOMT  by  iteratively  substituting  two  points 
in  D by  the  point  maXo{p,q)  until  the  origin  remains.  Again,  the  pair  of  points  p,q 
are  chosen  to  maximize  ||  maXo{p,q)  11- 

Let  Dz  denote  the  set  D at  the  point  in  the  construction  of  the  SOMT  in  which 
the  last  possible  pair  of  points  p,q  having  ||  maXo{p,q)  ||>  z have  been  joined,  and 
\etW,  = {peD,:\\p\\>z}. 

Lemma  6.3  During  the  construction  of  the  SOMT,  for  any  z > 0,  there  can  be  at 
most  one  point  in  IV^  on  or  above  the  horizontal  line  y = (0,z)  (also,  there  is  at  most 
one  point  in  on  or  below  the  horizontal  line  y = (0,  —z),  at  most  one  point  on 
or  to  the  right  of  the  vertical  line  x = (2,0)  and  at  most  one  on  or  to  the  left  of  the 
vertical  line  x = (—2,0)). 

Proof: 

Let  p,q  £ W;c  with  j/p,y,  > z.  Let  j/p  > j/,.  Then  maXo{p,q)  = (x,yq)  for  some  value 
X.  Hence  | | (x,j/,)  | |>  2,  contradicting  the  fact  that  p,q  ^ W^.  A similar  proof 
works  for  the  other  3 lines.  □ 

The  following  lemma  is  an  adapted  version  of  lemma  4 in  [32]. 

Lemma  6.4  Let  pi  = (x,,  j/,)  be  a point  in  a SOMT  generated  by  the  heuristic.  If 
X,-  7^  0 and  yi  ^ 0,  then  there  must  exist  a sink  on  the  vertical  path  x = x,  (horizontal 
path  y = y,)  starting  at  pi . 

Proof: 

Assume,  w.l.g.,  that  p,  lies  in  the  first  quadrant.  The  proof  is  by  induction  on  the 
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number  of  descendants  of  p,  . If  p,  does  not  have  descendants  then  it  is  a sink,  and  the 
result  immediately  follows.  Suppose,  on  the  contrary,  that  p,  = maXo{pj,Pk)-  Ob- 
serve that  these  two  points  belong  also  to  the  first  quadrant.  Thus  x,-  = min{xj,Xk), 
say  Xk,  and  p*  > y,.  By  induction  there  is  a sink  on  the  vertical  path  x — Xk  start- 
ing at  pj.  The  fact  that  yt  > y,  implies  that  this  sink  is  also  on  the  vertical  path 
X — Xi  — Xk  staring  at  p,.  A similar  argument  can  be  made  for  the  y-coordinate.  □ 
Now  assume  that  the  points  p and  q lie  in  the  first  quadrant.  The  following  two 
lemmas  are  taken  directly  from  [32].  The  proofs  can  be  found  in  their  paper. 

Lemma  6.5  During  the  construction  of  the  SOMT , for  any  z > 0 and  p,q  £ W^, 

1.  Xp  X,  and  pp  ^ y„ 

2.  Xp  < X,  if  and  only  if  yp  > y,, 

3.  if  Xp  < Xq,  then  Xp  -|-  y,  < 2. 

Lemma  6.6  Let  p^q  ^ with  Xp  < x,.  Then  there  exists  a sink  either  on  the 
horizontal  path  between  p and  (xp,yp)  or  on  the  vertical  path  between  (x,,pp)  and  q. 

It  is  straight  forward  to  prove  similar  results  when  both  of  the  points  p and  q lie 
in  any  of  the  other  three  quadrants.  The  following  theorem  provides  an  upper  bound 
for  the  traffic  generated  for  the  multicast  tree  by  the  heuristic.  Its  proof  is  similar  to 
the  proof  of  the  upper  bound  provided  in  [32]. 

Theorem  6.1  1{S0MT)  < 21{RSM A). 
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Proof: 

First,  note  that  for  any  RSA  A,  1{A)  = FI  LJ.  We  know  that  \ n A |>| 

MC{D,z)  I for  any  RSA  A.  Thus  it  is  sufficient  to  show  that  | SOMT  D |< 
2 I MC{D,z)  I for  any  z. 

For  any  given  z,  order  the  points  in  into  pi,p2, . . . ,Pm(z)  circularly  as  follows: 
first  take  the  nodes  in  the  first  quadrant  and  order  them  by  increasing  values  of  the 
x-coordinates  (with  Xi  >0),  then  the  nodes  in  the  second  quadrant  by  decreasing 
order  of  the  x-coordinate,  and  so  on. 

It  can  be  shown  that  MC{D,  z)  must  contain  distinct  points  for  each  pair  of  points 
P2t-i,P2i,  for  i = 1,2, . . . , [m(z)/2j.  Consider  the  two  pairs  of  points  (p2,_i,p2,), 
and  (p2i+i,P2i+2),  for  any  given  i in  the  above  range.  Assume,  for  simplicity,  that 
these  points  lie  in  the  first  quadrant.  By  Lemma  6.6  there  must  exist  a sink  on 
the  horizontal  path  y = y2»-i  between  P2i-i  and  (x2,, y2(-i)  or  on  the  vertical  path 
X = X2t  between  and  p2,.  Hence  MC{D,  z)  must  contain  a point  whose  x- 

coordinate  lies  between  z — y2,_i  and  X2j.  Also,  there  is  a sink  on  the  horizontal  path 
y = y2i+i  between  p2,+i  and  (a;2i+2,  y2t+i)  or  on  the  vertical  path  x = X2,+2  between 
(a^2i+2,y2.+i)  and  p2,+2-  Hence  MC{D,z)  must  contain  a point  whose  x-coordinate 
lies  between  z — y2i+i  and  X2,+2-  Since  X2,  + y2i+i  < by  Lemma  6.5,  then  the  above 
two  points  are  different.  A similar  result  holds  no  matter  where  the  points  lie. 

Thus  contains  at  most  2 | [MC{D,z)\  | points.  This  implies  that  if  m{z)  is 
even,  then  the  SOMT  generated  by  the  heuristic  intersects  at  most  2 | MC{ D,  z)  \ 
times.  If  m(z)  is  odd,  then  consider  p\. 
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Case  1:  Xp,  ^ 0. 

By  Lemma  6.4,  there  is  a sink  s on  the  vertical  path  x = Xp,  starting  at  pi.  Thus 
MC{D,z)  must  contain  a point  whose  x-coordinate  is  between  0 and  Xp, . For  the 
rest  of  the  \ \ —I  points  in  there  correspond  at  least  (|  \ — 1)/2  distinct 

points  in  MC{D,  z),  and  it  is  easy  to  see  that  these  points  must  be  different  from  the 
point  in  MC{D,z)  corresponding  to  pi.  Thus  MC{D,z)  contains  at  least  | \ /2 

distinct  points. 

Case  2:  Xp,  = 0.  We  can  assume  now  that  contains  points  on  or  above  the 
points  (0,  z),  {z,  0),  (0,  — ar),  ( — 2, 0)  in  the  same  vertical  or  horizontal  lines.  Otherwise, 
we  could  order  the  points  in  starting  with  the  points  in  some  other  quadrant  and 
reduce  the  problem  to  Case  1. 

Assume,  w.l.g.,  that  there  are  2c  points  in  in  the  first  quadrant,  for  some 
positive  integer  c.  Note  that  since  now  we  are  assuming  that  | | is  odd,  at  least 

one  quadrant  must  have  an  even  number  of  points  in  W,.  Thus  MC{D,z)  must 
contain  at  least  c points  in  the  first  quadrant. 

Case  2a.  MC{D,  2)  contains  at  least  c T 1 points  in  the  first  quadrant. 

MC{D,z)  must  contain  at  least  [(|  | — 2c)/2j  distinct  points  in  the  other 

three  quadrants.  Thus  MC{D,z)  contains  at  least  [|  \ /2J  — c + c+  1 >|  | /2 

distinct  points. 

Case  2b:  MC{D,z)  contains  exactly  c points  in  the  first  quadrant. 

We  can  assume  here  that  ||  ||<  z.  Otherwise,  we  can  order  the  nodes  circularly  the  other 

way  around,  starting  with  Pm{2)  and  finishing  with  pi. 
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We  will  show  that  at  leaist  one  of  the  points  (0,  z)  or  (z,  0)  cannot  be  one  of  these 
c points  in  MC{D,z).  Consider  the  first  2c  — 2 points  in  W^  in  the  first  quadrant, 
not  including  pi.  There  must  be  c — 1 distinct  points  in  MC{D,z)  to  cover  them, 
and  none  of  them  is  (0,  z)  or  (z,  0).  Thus  one  of  these  two  points  is  not  in  MC{D,  z). 
Assume,  w.l.g.,  that  pi  ^ MC{D,  z).  Thus  between  and  0 there  is  a different 

point  in  MC{D,z)  required  to  cover  the  sink  on  or  above  Pm{z)-  For  the  rest  of  the 
I I —2c  — 1 points  in  there  must  be  (|  | —2c  — l)/2  distinct  points  in 

MC{D,z).  Thus  MC{D,z)  contains  at  least  [|  \ /2J  — c + c -f-  1 >|  \ /2 

distinct  points. 

Thus  contains  at  most  2 | MC{D,z)  | points,  and  therefore  the  SOMT 
generated  intersects  Lz  at  most  2 | MC{D,  z)  | times.  Since  z is  arbitrary  this  holds 
for  any  z,  and  the  theorem  is  proved.  □ 

The  implementation  presented  in  [32]  for  their  algorithm  can  be  adapted  for  our 
algorithm.  The  key  idea  is  that  Steiner  points  can  only  be  introduced  by  the  heuristic 
from  adjacent  points  in  W^.  To  implement  the  algorithm  a priority  queue  is  used 
containing  all  the  sinks  below  the  lines  in  and  all  the  j Wz  | — 1 potential  new 
Steiner  points  maXo{p^q)  for  each  pair  of  adjacent  points  p,q  E See  [32]  for 

details  of  the  implementation. 

To  evaluate  the  performance  of  the  algorithm  we  ran  simulation  programs  for  an 
8*8  mesh  both  for  the  optimal  and  the  heuristic  algorithm.  The  optimal  solution  for 
100  randomly  generated  groups  was  found  by  exhaustive  search  of  all  shortest-path 
trees,  and  compared  to  the  suboptimal  solution  found  by  the  heuristic  algorithm  for 


Traffic 


66 


the  same  multiccist  groups.  The  traffic  is  the  number  of  links  in  the  multiccist  tree. 
The  results  are  shown  in  the  following  graph.  To  get  the  worst  case  line  we  added 
to  the  optimal  traffic  for  each  number  of  destinations  the  highest  difference  found 
between  the  optimal  traffic  and  the  traffic  computed  by  the  heuristic.  It  can  be  seen 
that  even  in  the  worst  case  the  values  computed  are  much  closer  to  the  optimal  values 
than  to  the  upper  bound. 


Figure  6.3.  Optimal  Algorithm  vs  Heuristic  (N=64) 
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6.2  Multiple  Groups 

The  solution  of  this  problem  is  very  similar  to  the  solution  provided  for  the  single 
group  case.  We  will  give  a generalization  of  the  algorithm  to  find  the  optimal  root 
given  the  m*n  mesh  and  the  p groups.  Once  we  have  computed  the  optimal  root, 
the  heuristic  algorithm  proposed  in  the  previous  section  can  be  used  to  construct  the 
SOMT. 

Now  we  want  to  find  the  node  i such  that 

p 

c,  = ^ max  d[i,Tijk) 

, 1 <!  A;  <C  m , 

J 

is  minimized,  where  rijk  is  the  kth  element  of  the  jth  destination  list,  and  mj  is  the 
number  of  nodes  in  the  jth  group.  Let  maxi<t<„i^  d{a,  rijk)  = d{a,  ruja),  for  any  node 
a.  nrija  is  now  the  node  in  group  j farther  away  from  node  a. 

Lemma  6.7  Let  y be  an  optimal  solution.  For  any  node  x,  if  Cy  < c^,  then  x has  a 
neighbor  or  a diagonal  z such  that  < c^. 

Proof: 

As  before,  let  x = (f,;),  y = {k,l).  Assume,  w.l.g.,  that  x and  y differ  in  their  first 
component,  with  i < k,  and  that  j < 1.  Again,  define  Zi  = {i  + 1,;),  = {i,j  + 1), 

and  Z3  = (f  + 1,  J + 1).  Let  be  the  cost  of  node  n when  only  group  j is  considered 
(as  a single  group). 

Case  1:  < Cjj,  for  some  j,  I < j < p.  An  immediate  corollary  of  Lemma  6.2 

is  that  Cjz3  < Cjr  VI  < j < p.  Thus  < c^,  and  the  lemma  is  true. 
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Case  2:  Vj,  I < j < p-  Thus,  < Cjx  for  i=l  or  2.  Note  also  that 

I ^jz,  ~ Cjx  |<  1. 

If  p is  odd,  then  < Cjx  for  at  least  fp/2]  groups  for  either  i=l  or  i=2.  Thus 
Cz,  < Cx  for  this  value  of  i,  and  we  are  done.  If  p is  even,  we  may  think  that  maybe 

^ Cjx  for  p ! groups,  but  for  the  other  p^2  groups,  Cj^^  ^ Cjx^  and  that  the  same 
thing  occurs  with  ^2,  and  thus  none  has  a cost  lower  than  x.  We  will  show  that  this 
cannot  happen.  We  just  need  to  show  that  for  2 groups,  the  Ccise  in  which  < Cix, 
Ciz2  > cix,  and  C2z,  > C2x,  C2z2  < C2r  cannot  occur. 

Assume  that  the  above  4 inequalities  hold.  The  first  2 inequalities  imply  that  mix 
has  first  component  > i+1.  It  must  have  second  component  < j,  since  = cix-  The 
other  2 inequalities  imply  that  m2i  heis  first  component  < i,  and  second  component 
> J + 1,  as  shown  in  figure  6.4. 

But  we  are  assuming  that  = Cjx  for  j=l,2.  Thus  d{y,mjx)  < d{zz,mjx)  for 
j = l,2.  This  is  impossible  given  the  locations  of  mix  and  m2x  in  the  figure.  Thus  we 
have  a contradiction,  and  thus  it  must  be  the  case  that  zi  or  22  have  lower  cost  than 

X.  □ 

The  algorithm  to  find  the  root  of  the  SOMT  is  thus  essentially  the  same  as  for  the 
case  of  a single  group.  First  we  find  the  node  that  is  in  the  center  of  the  smallest  grid 
that  contains  all  the  nodes  in  all  the  destination  lists,  and  from  there  if  necessary  we 
move  to  lower  cost  nodes  until  we  converge  to  the  optimal  root. 

In  the  single  group  case,  the  root  of  the  tree  is  always  the  primary  destination 
of  the  group.  When  we  have  multiple  groups,  we  need  to  compute  the  primary 
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destination  for  each  group.  Now  a way  to  compute  the  primary  destinations  while 
the  routing  is  being  determined  is  presented. 

6.3  Primary  Destinations 

The  following  algorithm  computes  the  primary  destinations  of  the  p groups  while 
determining  the  routing  using  the  heuristic  algorithm  presented  previously. 

For  a given  point  p generated  by  the  algorithm,  ccp  represents  the  set  of  multiccist 
groups  containing  a member  of  the  group  that  is  a descendent  of  p in  the  SOMT.  For 
each  group  i,  NC{  denotes  the  number  of  disjoint  subtrees  containing  a member  group 
i.  When  two  points  p,q  are  joined  and  substituted  by  maXo{p,q),  every  descendent 
of  p and  q becomes  a descendent  of  maXo{p,q).  Let  p,  denote  the  nodes  in  the  I'th 
group.  Initially,  cCp  includes  i if  p G p,.  Also,  if  a group  i does  not  have  any  member 
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Figure  6.4. 
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that  is  a descendent  of  p or  q,  but  maXo{p,  q)  € gi,  then  i must  be  added  to  cCmaxo(p,q)- 
When  all  the  nodes  in  group  i belong  to  the  tree,  and  NCi  = 1,  implying  that  all  of 
them  are  descendants  of  the  same  point,  then  the  primary  destination  for  this  group 
can  be  computed.  It  must  be  the  Icist  node  included  in  the  tree,  i.e.,  maXo{p,q).  Let 
D now  represent  the  set  of  distinct  nodes  in  the  multicast  groups. 


Algorithm  6.2  Primary  destinations  for  mesh  networks 
begin 

cCp  = <l>  for  every  node  p; 
ccp  = CCp  U i'll  p £ p,; 

compute  optimal  root; 
while  I Z)  |>  1 do 

let  p,q  be  substituted  by  maXo{p,q)] 

^^maxo(p,q)  — bl  CCg, 

V group  i not  marked  do 
if  maXo{p,  q)  G p,  then 

mark  maXo{p^  q)  in  p,  as  visited; 

('^^axo(p,g)  ~ ^^maxo(p,q) 

if  NCi  = 0 and  maXo{p,  q)  G p,  then  NCi  = 1 
else  if  i G cCp  and  i € cc,  then  NCi  = NCi  ~ 1 

else  if  i ^ cCp  and  i ^ cc,  and  i G cCjnaxo{p,q)  then  NCi  = NCi  -f  1; 
if  NCi  = I and  all  nodes  in  group  i are  marked  as  visited  then 
PDi  = maXo{p,q)-, 
mark  group  i; 

end; 


6.4  Multiple  Trees 

First  the  subtrees  for  the  connected  components  of  graph  G'  as  defined  previously 
must  be  built.  These  subtrees  contains  the  overlapping  nodes,  and  by  maintaining 
the  relative  position  of  the  nodes  in  each  subtree  the  multiple  group  ordering  property 


is  satisfied. 
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Now  we  examine  how  to  construct  the  multicast  tree  T,  for  group  i.  Suppose  that 
the  subtrees  that  must  be  part  of  T,  have  been  determined.  As  said  before,  a depth 
equal  to  the  maximum  distance  from  the  root  of  each  subtree  to  the  farthest  node 
in  the  subtree  that  is  a member  of  group  i is  assigned  to  each  of  these  subroots.  It 
may  be  possible  that  the  same  node  is  root  of  two  or  more  subtrees.  In  that  case 
the  maximum  of  the  depths  is  assign  to  it.  The  nodes  in  group  i that  did  not  belong 
to  any  of  the  connected  components  but  that  can  be  inserted  in  any  of  the  subtrees 
without  increasing  the  cost  of  the  subtree  are  included  in  the  subtree.  These  nodes 
could  increase  the  traffic  generated  by  the  subtree,  but  not  its  cost.  The  rest  of  the 
nodes  in  group  i are  assigned  a depth  of  0.  Now  these  nodes  need  to  be  connected 
together  and  to  the  roots  of  the  subtrees.  We  have  a single  group  with  a depth 
assigned  to  each  node,  and  an  optimal  root  must  be  found.  A similar  algorithm  as 
the  one  used  to  find  the  optimal  root  for  a mesh  network  when  a single  group  was 
given  can  be  used.  The  cost  of  node  i is  now  given  by  c,  = max j^o {d{i,  j) -[■  dept h{j)} . 
The  same  proof  discussed  in  section  6.1  works  here.  If  a node  is  not  an  optimal  root, 
then  it  must  have  a neighbor  with  lower  cost.  Thus  the  same  algorithm  used  before 
to  converge  to  an  optimal  root  moving  in  each  iteration  to  the  direction  of  a lower 
cost  node  can  be  used  as  well. 

Simulations  ran  for  an  8*8  mesh  network  under  the  same  conditions  as  the  ones 
ran  for  linear  and  ring  networks  show  similar  results.  The  multiple  trees  approach 
shows  to  be  slightly  better  in  terms  of  the  cost  of  the  multicast  and  much  better  in 
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terms  of  the  variability  of  the  primary  destinations,  while  it  generates  more  traffic 
than  the  single  global  tree  approach. 


Cost  of  Multicast 
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Figure  6.5.  Cost  for  Mesh  Network  (N=64) 


Traffic 
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Figure  6.6.  Traffic  for  Mesh  Network  (N=64) 


% of  Identical  Primary  Destinations 
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Figure  6.7.  Primary  Destinations  for  Mesh  Network  (N=64) 


CHAPTER  7 

MULTICAST  TREES  FOR  HYPERCUBE  NETWORKS 
7.1  Single  Group 

Let  Qn  be  a hypercube  multiprocessor  with  N = 2'^  nodes,  and  let  D he  & set 
of  points  from  the  hypercube.  A node  x is  labeled  by  the  binary  number  b[x)  = 
6„_i6„_2  . . . 6q.  The  distance  d in  Qn  between  two  nodes  is  equal  to  the  Hamming 
distance  of  their  binary  addresses. 

The  first  step  to  construct  the  multicast  tree  for  the  set  D is  to  find  the  optimal 
root,  that  is,  the  node  i that  minimizes  c,  = maxj^£)d(i,j). 

Lemma  7.1  Let  Q,  be  the  smallest  subcube  of  Qn  that  contains  all  the  nodes  in  D, 
and  let  x be  an  optimal  root.  Then  x ^ Qs- 

Proof: 

Assume,  w.l.g..,  that  the  first  n — s bits  of  the  nodes  in  D are  all  the  same,  say  all 
ones.  Assume,  contrary  to  the  statement  of  the  Lemma,  that  x is  an  optimal  root 
but  X ^ Q^.  Then  at  least  one  of  the  first  n — s bits  of  x is  0.  Let  i be  one  of  those 
bit  positions.  Consider  node  y with  all  bits  cis  x except  bit  position  z,  in  which  y has 
a 1 instead  of  a 0.  Clearly,  d[y,a)  < d[x,a)  Va  € D.  Thus  Cy  < c^,  contradicting  the 
fact  that  X is  optimal.  Thus  x 6 Qs,  and  the  Lemma  is  proved.  □ 
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Let  D = {ni,  U2, . . . , n^},  where  m =\  D \,  and  let  be  the  smallest  subcube 
of  Qn  containing  all  the  nodes  in  D.  Let  n7  denote  the  complement  of  node  n,,  i.e., 
the  node  obtained  by  reversing  every  bit  in  the  binary  address  of  n,  not  including 
the  bit  positions  in  which  all  the  nodes  in  D agree  upon.  Note  that  to  form  n~  only  s 
bits  are  complemented,  and  that  d(n,,  u7)  = s.  The  cost  of  the  complement  of  every 
node  in  D is  s.  The  cost  of  the  neighbors  of  each  complement  of  every  node  in  D 
is  at  least  s — 1,  and  so  on.  This  is  the  idea  behind  the  following  algorithm,  which 
computes  the  optimal  root  in  0{N' logN')^  where  N'  = 2". 

Algorithm  7.1  Optimal  Root 

begin 

p={}: 

for  i=0  to  — 1 do  c,  = — 1; 
for  i=l  to  m do 
Cn7  — 

p = pn  {n-}; 

while  there  are  nodes  with  cost  -1  do 
q = p; 
p=  {}; 

Vy  € 9 do 

V neighbor  i of  y do 
if  Cx  = —1  then 

Cx  — Cy  1; 

p = pn  {i}; 

every  node  in  p is  an  optimal  root 
end; 

The  complexity  of  this  algorithm  may  be  too  expensive,  especially  if  there  are 
only  a few  nodes  in  the  multicast  group  and  s = n.  We  propose  a heuristic  algorithm 
to  find  a suboptimal  root  more  efficiently.  The  main  result  used  to  find  the  optimal 
root  for  mesh  networks  does  not  always  produce  the  optimal  root  for  hypercubes.  A 
node  which  is  not  optimal  may  have  all  its  neighbors  with  higher  cost  than  itself. 
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However,  the  same  idea  can  be  used  as  a heuristic  to  find  a suboptimal  root.  The 
heuristic  provides  a better  result  if  we  start  with  a good  candidate  for  optimal  root. 

Let  n,j  denote  the  jth  bit  of  node  n,.  Define  node  x in  the  following  way:  Xj  = 1 
iff  Tiij  > m/2.  Node  x hcis  a one  in  bit  position  j if  and  only  if  more  than 
half  of  the  nodes  in  the  group  have  ones  in  that  bit  position.  This  node  is  called 
the  median  of  D.  In  average,  we  should  expect  the  optimal  root  to  be  close  to  the 
median.  Thus,  the  median  x of  £)  is  our  initial  candidate  for  the  root.  Then  the 
neighbors  of  x are  examined  to  see  if  any  of  them  has  a lower  cost  than  x.  In  this  case 
X is  replaced  by  this  neighbor  and  the  process  continues.  Eventually  the  algorithm 
converges  to  a node  that  does  not  have  a neighbor  with  lower  cost.  The  cost  of  the 
median  can  be  computed  in  0{mn).  Once  the  cost  of  a node  is  computed,  the  cost 
of  the  neighbors  can  be  computed  in  0(m).  At  most  n hops  are  required  to  converge 
to  the  final  root,  so  the  total  complexity  of  the  algorithm  is  0{mn). 

We  ran  simulations  for  hypercubes  going  from  2 to  2048  nodes  for  randomly 
generated  multicast  groups.  Each  node  had  a fixed  probability  of  0.2  of  being  included 
in  the  multicast  group.  The  table  below  shows  the  results  obtained.  N denotes  the 
number  of  nodes  in  the  hypercube.  The  second  component  shows  the  percentage  of 
CcLses  in  which  the  root  computed  by  the  algorithm  is  in  fact  the  optimal  root.  The 
low  average  in  the  number  of  hops  required  showed  in  the  third  column  shows  that 
in  general,  for  hypercubes  of  these  sizes,  the  algorithm  converges  to  the  root  much 
faster  than  the  O(mn^)  worst  case  time  complexity.  It  is  also  important  to  mention 
that  in  all  the  cases  (except  for  one  of  the  runs  for  N = 64)  that  the  heuristic  did 
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not  compute  the  optimal  root,  the  difference  between  the  cost  of  the  optimal  and  the 
cost  of  the  node  computed  by  the  algorithm  was  only  one.  Thus  it  is  fair  to  say  that 
the  heuristic  provides  a very  good  approximation  for  the  optimal  root. 


N 

% Optimal  Root 

Avg  of  Hops 

2 

100 

0 

4 

100 

.04 

8 

100 

.11 

16 

99.2 

.22 

32 

96.9 

.42 

64 

91.9 

.53 

128 

90.5 

.62 

256 

90.9 

.61 

512 

93.8 

.79 

1024 

89.0 

.76 

Choi  and  Esfahanian  [6]  have  shown  that  the  OMT  problem  is  NP-Complete  for 
hypercubes  networks.  Lan  et  al.  [28]  proposed  a heuristic  algorithm  {LEN  heuristic) 
to  construct  a sub-optimal  multicast  tree  for  hypercubes.  In  their  heuristic,  the  root 
of  the  tree  is  the  source  of  the  multicast  message. 

Each  node  i in  the  tree  (beginning  with  the  source)  executes  a multicast  algorithm 
to  select  the  adjacent  nodes  to  which  it  will  send  the  message.  The  multicast  group 
is  partitioned  into  disjoint  multicast  subgroups,  each  consisting  of  destination  nodes 
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which  will  be  descendants  of  node  i in  the  tree.  The  multiccist  algorithm  selects 
an  adjacent  node  with  a maximum  number  of  descendants  in  the  multiccist  group  it 
received. 

Consider  the  hypercube  with  2"  nodes.  Assume  each  node  v is  labeled  by  the 
binary  number  b{v)  = 6n_i6„_2  ...  to-  The  relative  address  of  a node  t(u,)  with  respect 
to  a node  b{vj),  denoted  by  6(u,ju_,),  is  defined  by  taking  the  bit-wise  exclusive-or  of 
the  two  nodes.  Let  X be  a collection  of  m n-bit  binary  numbers  a*_j,  a^_2, . . . , ao> 
1 < A:  < m,  and  let  B{X)  = X„_2, . . . , A'o),  where  X{  — 5Zr=i  The 

multica.st  algorithm  executed  by  node  v upon  receiving  the  multicast  subgroup  X,  is 
cis  follows: 

1.  For  each  node  b{x)  6 X,  replace  b{x)  by  b{x\v). 

2.  Compute  B{X)  = . . . , A^o),  and  select  the  smallest  i such  that  A,  > Xj, 

for  any  j. 

3.  If  A,  = 0,  stop.  Otherwise,  create  a multicast  subgroup  Y containing  all 
b{x)  € A such  that  b{x)  has  a one  in  bit  position  i. 

4.  Send  Y to  the  node  in  the  hypercube  adjacent  to  v such  that  6(u|u)  has  a one 
in  bit  position  i. 

5.  Set  X — X — Y , and  go  to  step  2. 

A problem  with  this  algorithm  is  that  it  often  splits  neighboring  nodes  among 
different  multicast  subgroups,  which  increases  the  traffic  generated.  Choi  and  Es- 
fahanian  have  proposed  a modification  of  this  algorithm  (which  they  call  Covered 
heuristic)  that  keeps  pairs  of  adjacent  nodes  together  in  the  same  subgroups. 
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A node  x in  a multicast  group  or  destination  sublist  is  covered  by  an  adjacent 
node  y in  Qn  if  d(x,r)  > d(y,r),  where  r is  the  root  of  the  multicast  tree  being 
generated.  The  Covered  heuristic  is  as  follows: 

1.  For  each  node  6(x)  6 X,  replace  b(x)  by  fc(x|u).  Also,  mark  each  node  in  X as 
covered  or  not  covered. 

2.  Compute  B{X)  = (X„_i , . . . , Xq).  Let  / be  the  set  of  integers  z,  such  that 
X{  > Xj  \/j  = 0,1, ...,n  — 1.  If  I 7 1=  1 then  go  to  step  3.  Otherwise,  let  W = 
{b{x)  G X I b{x)  is  not  covered}.  Compute  B{W),  and  select  the  smallest  i such  that 
Wi  > Wj  Vi  = 0,l,...,n-1. 

3.  If  Xi  = 0,  stop.  Otherwise,  create  a multicast  subgroup  Y containing  all 
b{x)  G X such  that  b{x)  has  a one  in  bit  position  i. 

4.  Send  Y to  the  node  in  the  hypercube  adjacent  to  v such  that  fc(u|i;)  hcis  a one 
in  bit  position  i. 

5.  Set  X = X — Y , and  go  to  step  2. 

The  changes  in  the  second  step  of  the  algorithm  attempt  to  keep  pairs  of  adjacent 
nodes  together  in  the  same  sublists. 

Example  7.1  : 

Consider  the  following  multicast  group:  {0100,0101,1010,1011}.  An  optimal  root 
for  this  group  is  node  0000.  The  multicast  trees  generated  by  each  of  these  heuristics 


are: 
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0000 


0001  0010  0100 


0011  0101  1010 
1011 


0000 

/\ 

0100  0010 

0101  1010 
1011 


LEN  Heuristic 


Covered  Heuristic 


In  this  case  the  nodes  in  D 0101  and  1011  are  covered.  The  Covered  heuristic  kept 
these  two  nodes  together  with  the  nodes  that  covered  them,  while  the  LEN  heuristic 
split  the  pairs  of  nodes  0100,  0101,  and  1010,  1011.  However,  the  Covered  heuristic 
was  able  to  keep  these  two  pairs  of  nodes  together  in  the  same  sublists  only  because 
I / |>  1 in  the  second  step  of  the  algorithm.  It  may  happen  that  some  nodes  could 
be  covered  by  adjacent  nodes  in  the  destination  sublist,  but  they  get  split  because 
I / 1=  1.  Consider,  for  example,  adding  the  node  0001  to  the  multicast  group  in 
the  previous  example.  Now  B{X)  = {2, 2, 2, 3}  in  the  first  pass  of  the  algorithm,  so 
I / 1=  1,  and  the  following  tree  is  generated: 


0000 


0001  0010  0100 


0011  0101  1010 


1011 
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Nodes  1010  and  1011  are  now  split,  generating  a traffic  of  7 units.  Had  we  applied 
the  covered  idea  even  when  | / |=  1 the  following  tree  would  have  been  generated: 


0000 


0100  0010  0001 
0101  1010 

1011 


with  a traffic  of  only  6. 

This  motivates  the  following  modification  o the  Covered  heuristic:  in  the  second 
step  of  the  algorithm,  compute  B{W)  and  the  corresponding  value  of  i without 
computing  set  1.  The  new  algorithm  produces  the  optimal  multicast  tree  shown  in 
the  previous  figure  for  this  example. 

The  performance  of  this  algorithm  was  compared  to  the  Covered  heuristic  for 
randomly  generated  multicast  groups  for  hypercubes  going  from  64  nodes  to  1024 
nodes.  The  results  are  shown  in  the  following  graph. 


Traffic 
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Number  of  Nodes 


Figure  7.1.  Covered  heuristic  vs  Our  heuristic 
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7.2  Multiple  Groups 

The  solution  to  this  problem  is  similar  to  the  solution  provided  for  the  single 
group  case.  The  algorithm  to  find  an  optimal  root  given  Qn  and  a multicast  group 
D can  be  generalized  for  p groups.  When  a single  group  Wcis  given  the  cost  of  each 
node  Wcis  computed.  The  algorithm  can  be  applied  for  each  of  the  p groups,  and  the 
node  i that  minimizes  C{  = 5Zi=i  rnaxi<fc<m_,  d{i,njk)  can  be  found.  The  values  given 
by  maxi<jt<mj  d{i,njk)  are  the  values  computed  by  the  algorithm  for  each  node,  so 
the  node  that  minimizes  the  above  sum  can  be  found  with  Np  computations.  Thus 
the  optimal  root  for  p groups  has  0{NptogN)  time  complexity.  A similar  algorithm 
as  the  heuristic  algorithm  presented  before  can  also  be  used  for  multiple  groups  to 
find  a suboptimal  root. 

The  multicast  tree  for  the  p groups  can  be  generated  using  the  same  heuristic  used 
for  a single  group  which  constructs  a shortest-path  tree  and  attempts  to  minimize 
the  traffic  generated.  The  primary  destination  for  each  group  is  computed  while  the 
routing  is  being  determined. 

7.3  Primary  Destinations 

The  following  algorithm  computes  the  primary  destinations  of  the  p groups  while 
the  routing  for  the  multicast  tree  is  being  generated.  The  heuristic  algorithm  pre- 
sented in  section  7.1  is  used.  The  primary  destination  for  each  group  i is  initialize  to 
be  the  root  of  the  tree.  Each  time  a new  sublist  is  generated  by  the  algorithm,  the 
primary  destination  for  group  i can  be  set  to  be  the  receiver  of  the  sublist  (node  u 
in  the  algorithm)  if  and  only  if  the  sublist  is  a subset  of  p,.  If  the  sublist  contains 
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some  nodes  in  and  some  that  do  not  belong  to  then  PDi  cannot  change  any 
more.  Variable  donci  is  set  to  true,  meaning  that  this  group  does  not  have  to  be 
considered  any  more  for  the  computation  of  the  primary  destination.  The  algorithm 
is  as  follows: 

1.  For  each  group  i,  initialize  PD{  = root,  donci  = false 

2.  For  each  node  b{x)  € X,  replace  b{x)  by  6(i|u).  Also,  mark  each  node  in  X as 
covered  or  not  covered. 

3.  Let  W = {b{x)  € X \ b{x)  is  not  covered}.  Compute  B{W)  = {Wn-\,  Wn-2,  ■ ■ ■ , fko), 
and  select  the  smallest  i such  that  IF,  > Wj  Vj  = 0,l,...,n  - 1. 

4.  If  Xi  = 0,  stop.  Otherwise,  create  a multicast  subgroup  Y containing  all 
b{x)  G X such  that  6(i)  has  a one  in  bit  position  i. 

5.  For  each  group  with  donci  = false  do 

if  PDi  = V then 
if  ga  G Y then 

if  9i  Q Y then  PDi  = u (cis  in  step  6) 
else  donci  = true 

6.  Send  Y to  the  node  u in  the  hypercube  adjacent  to  v such  that  i(u|u)  has  a 
one  in  bit  position  i. 

1.  Set  X = X — Y , and  go  to  step  3. 

7.4  Multiple  Trees 


As  before,  the  subtrees  for  the  connected  components  of  graph  G'  as  defined  pre- 
viously must  be  built.  Once  the  subtrees  that  must  be  part  of  the  multicast  group  for 
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each  group  i have  been  determined,  a depth  equal  to  the  maximum  distance  from  the 
root  of  each  subtree  to  the  farthest  node  in  the  subtree  that  is  a member  of  group  i is 
assigned  to  each  of  these  subroots.  The  nodes  in  group  i that  did  not  belong  to  any 
of  the  connected  components  but  that  can  be  inserted  in  any  of  the  subtrees  without 
increasing  the  cost  of  the  subtree  are  included  in  the  subtree.  The  rest  of  the  nodes 
in  group  i are  cissigned  a depth  of  0.  Now  these  nodes  need  to  be  connected  together 
and  to  the  roots  of  the  subtrees.  We  have  a single  group  with  a depth  assigned  to 
each  node,  and  an  optimal  root  must  be  found.  The  following  algorithm  computes 
the  optimal  root  for  the  multicast  tree  of  group  i: 

Copt  = 0;  opt. root  = 0; 
for  i=0  to  N-1  do 
c,  = 0; 

for  j=l  to  m do 

if  c,  < d{i,j)  + depth{j)  then  c,  = d{i,j)  + depth(j)-, 
if  c,  < Copt  then 

^opt  — 

optjroot  = i] 

The  complexity  of  this  algorithm  is  0{mN),  where  m = Yl^=i  ^.nd  m,-  is  the 
number  of  nodes  in  group  i.  A more  efficient  algorithm  can  be  used  to  find  a subopti- 
mal  root  using  a similar  heuristic  cis  the  one  discussed  in  the  previous  section.  It  can 
be  modified  to  favor  nodes  with  higher  depth.  Let  5,  denote  the  set  of  nodes  in  set  A to 
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be  connected  and  depth{i)  the  corresponding  depths.  The  median  of  the  groups  is  now 
defined  ais  the  node  x where  x,-  = 1 f//  gi*{depth{i)  + l)  > Xllli(dept/i(i)  + l)/2. 

The  routing  is  generated  with  the  same  heuristics  as  discussed  in  the  section 
7.1.  Simulations  ran  for  Qq  under  the  same  conditions  ais  the  ones  ran  for  the  other 
network  topologies  discussed  in  previous  chapters  show  similar  results.  The  multiple 
trees  approach  shows  to  be  slightly  better  in  terms  of  the  cost  of  the  multicast  and 
much  better  in  terms  of  the  variability  of  the  primary  destinations,  while  it  generates 
more  traffic  than  the  single  global  tree  approach. 


Cost  of  Multicast 
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Figure  7.2.  Cost  for  Hypercube  Network  (N=64) 


Traffic 
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Number  of  Groups 


Figure  7.3.  Traffic  for  Hypercube  Network  (N=64) 


% of  Identical  Primary  Destinations 
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Number  of  Groups 


Figure  7.4.  Primary  Destinations  for  Hypercube  Network  (N=64) 


CHAPTER  8 

CONCLUSION  AND  FUTURE  WORK 


The  objective  of  our  work  was  to  provide  message  ordering  in  a multicast  envi- 
ronment given  a set  of  multicast  groups.  We  have  proposed  multicast  trees  to  achieve 
this  goal.  For  a single  multicast  group  a single  tree  is  a sufficient  and  efficient  way 
of  providing  message  ordering.  The  root  of  the  tree  can  provide  a total  order  for  all 
the  nodes  in  the  group.  When  multiple  groups  are  given,  two  graph  models  can  be 
used  to  provide  the  multiple  group  ordering  property:  a single  global  tree  can  be  con- 
structed for  all  the  nodes  in  overlapping  groups,  or  separate  trees  can  be  constructed 
for  each  group  with  common  subtrees  for  overlapping  nodes.  The  goal  in  both  cases 
is  to  build  trees  which  minimize  time  and  traffic. 

Algorithms  to  build  optimal  multicast  trees  have  been  discussed  for  linear  and  ring 
networks.  The  problem  of  building  optimal  trees  for  hypercubes  has  been  shown  to  be 
an  NP-Complete  problem.  For  mesh  networks  the  problem  is  conjectured  to  be  also 
NP-Complete.  For  these  topologies,  heuristic  algorithms  have  been  proposed  that 
build  shortest-path  trees  while  keeping  the  traffic  generated  low.  For  mesh  networks, 
it  remains  an  open  problem  to  determine  if  there  is  a polynomial  time  algorithm  to 
construct  an  OMT.  Some  similar  Steiner  tree  problems  have  been  shown  to  be  NP- 
Complete,  but  their  proofs  do  not  seem  to  be  applicable  to  the  OMT  problem.  The 
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rectilinear  Steiner  tree  problem  is  also  an  important  problem  in  determining  global 
routing  and  wiring  for  VLSI  circuit  layout. 

For  all  the  network  topologies  discussed  in  this  work,  the  comparison  made  be- 
tween the  single  tree  approach  and  the  multiple  trees  approach  revealed  similar  re- 
sults. The  single  global  tree  for  all  overlapping  groups  showed  a slightly  higher  cost 
for  the  multicasts  but  generated  less  traffic  than  the  multiple  trees  technique.  How- 
ever, the  bigger  difference,  as  expected,  was  in  the  primary  destinations.  The  multiple 
trees  approach  showed  much  more  variability  in  the  primary  destinations. 

Some  related  problems  were  not  addressed  in  this  work.  For  example,  the  im- 
pact that  different  switching  strategies  as  store-and-forward,  virtual  cut-through  and 
wormhole  may  have  on  the  performance  of  the  system  when  used  together  with  the 
routing  strategies  was  not  considered.  Perhaps  the  major  weakness  of  our  approach 
to  provide  message  ordering  is  the  cost  of  building  the  multicast  trees.  They  should 
be  used  for  multicast  groups  whose  life  is  relatively  long,  so  the  initial  cost  of  build- 
ing the  trees  is  distributed  among  many  messages.  Only  static  multicast  groups  were 
considered  in  this  work.  Methods  to  modify  the  algorithms  to  allow  for  deleting, 
adding  or  modifying  groups  dynamically  should  be  further  studied.  Another  future 
directions  of  this  work  is  to  consider  how  to  incorporate  reliability  issues  into  our 
scheme.  All  these  problems  need  to  be  addressed  to  develop  a complete  scheme  to 


be  implemented  in  a real  network. 
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